NASA-CR— 19516^ 



c> •- ■ x 


mP 


ELECTROMAGNETICS LABORATORY 
TECHNICAL REPORT 94-4 

October 1993 


INTEGRATED REFLECTOR ANTENNA DESIGN AND ANALYSIS 


FINAL REPORT 


[. L. Zimmerman, S. W. Lee, S. Ni, M. Christensen and Y. M. Wang 
University of Illinois, Urbana-Champaign ® 




pTri/l 




Supported by 

NASA Lewis Research Center 
Cleveland, Ohio 

NCC3-156 


Electromagnetics Laboratory 
Department of Electrical and Computer Engineering 
University of Illinois at Urbana-Champaign 
Urbana, Illinois 61801 


o 

c 

o o — 

UJ Z- 
< — 

< *-• 
a Z ^ 

O O 
UJ ' 
y— i/) 

T uj u 

H O O 

a 

^ i” O' 

z 

NO UJ 


0^ < •- ro 
u ro 

I * 

O' O co 

O ^ 

| U ^ ^ 
<t UJ > • 

oO -J _J > 
<f U_ < — 

U> z c 
o 


PW&CSCHN6 PAGE W ANK NOT FILMED 


u 


INTEGRATED REFLECTOR ANTENNA DESIGN AND ANALYSIS 


M. L. Zimmerman, S. W. Lee, S. Ni, M. Christensen and Y. M. Wang 
University of Illinois, Urbana-Champaign 


Final Report 

NASA NCC3-156 


October 1993 


Electromagnetics Laboratory 
Department of Electrical and Computer Engineering 
University of Illinois at Urbana-Champaign 
Urbana, IL 61801 



ABSTRACT 


Reflector antenna design is a mature field and most aspects have been studied. 
However, of that most previous work is distinguished by die fact that it is narrow in scope, 
analyzing only a particular problem under certain conditions. Methods of analysis of this 
type arc not useful for working on real-life problems since they can not handle the many 
and various types of perturbations of basic antenna design. In this thesis, the idea of an 
integrated design and analysis is proposed. By broadening the scope of the analysis, it 
becomes possible to deal with the intricacies attendant with modem reflector antenna design 
problems. 

In this diesis, the concept of integrated reflector antenna design is put forward. A 
number of electromagnetic problems related to reflector antenna design arc investigated. 
Some of these show how tools for reflector antenna design arc created. In particular, a 
mctiiod for estimating spillover loss for open-ended waveguide feeds is examined. The 
problem of calculating and optimizing beam efficiency (an important figure of ment in 
radiomciry applications) is also solved. Other chapters in this thesis deal with applications 
of this general analysis. The wide-angle scan abilities of reflector antennas is examined and 
a design is proposed for the ATDRSS triband reflector antenna. The-following chapter 
discusses the development of a general phased-array pattern computation program and 
shows how the concept of integrated design can be extended to other types of antennas. 

The conclusions are contained in the final chapter. 
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1. INTRODUCTION 

The subject of this thesis is based on the analysis and design of reflector antennas 
using an integrated approach. The concept is to obtain a method of analysis that can be 
applied to a wide variety of problems, rather than one that is narrow in scope and may be 
applied to only a few specific problem types. Reflector antenna design is a mature field. 
Reflector antennas have been in common use since World War EL They are still popular 
today due to their simplicity, ease of construction, low cost, and light weight Most 
aspects of reflector antenna design have been studied in the past Figures of merit such as 
directivity, beam efficiency, and sidelobe level have been looked at closely. The effect of 
feeds, compensation for reflector distortion, shaping of the reflector dish, etc. have also 
been studied. What most of these earlier studies lacked, however, was a broad scope. A 
partic ular aspect of reflector antenna design was studied. A method of analysis was 
derived and results were obtained. Then another problem would be tackled, with its own 
solution. However the method of analysis could rarely be used for anything other than 
generating research papers on that particular topic. 

The intention of this thesis is to make the jump from methods of analysis that 
generate research papers to methods of analysis that can be used to tackle a broad range of 
real life design problems. This means that the anaysis must be able to handle a broad range 
of reflector antenna intricacies. The method of analysis must be very general in nature and 
yet be simple enough so that it may be used by someone other than the originator of the 
method of analysis. 

In each chapter a separate topic of reflector antennas is discussed. The chapters are 
related by the fact that they represent some of the problems encountered with reflector 
design. The second chapter discusses an important parameter for reflector antenna design, 
the spillover loss. In the third chapter, the aperture integration method, which is used in 
our reflector antenna analyses, is explained in detail. The fifth chapter examines the 
analysis of beam efficiency, the most important figure of merit for antennas used in 



radiometry applications. The fourth and sixth chapters look at problems posed for reflector 
antennas. In the fourth chapter, a study is undertaken to find the limits to which a reflector 
an t*™™ system might be used for scanning. In the sixth chapter the challenge of 
constructing a triband reflector antenna is examined. In this chapter, in particular, all the 
capabilities of the method of analysis must be brought to bear on the problem, since there 
are many constraints upon the design. It is this type of reflector a n ten n a problem, 
involving a space-borne antenna, that will be the focus of much of the future research in 
reflector antennas. Conclusions are presented in the final chapter. 



2 REFLECTOR SPILLOVER LOSS OF AN OPEN-ENDED RECTANGULAR AND 

CIRCULAR WAVEGUIDE FEED 


Open-ended rectangular and circular wavegides are a most commonly used feed for 
reflector antennas 1 . Their radiation pattern was first calculated by Chu in 1940 by using a 
form of KirchhofF s approximation [1,2]. Important pattern characteristics such as 3 dB 
beamwidth, zeros, and sidelobes are well documented in the literature [3,4]. In reflector 
applications, there is another important parameter, namely, the spillover loss, which is 
as the fraction of power received by the symmetrical reflector within the half-cone 
angle 0 (Figure 2.1). (All figures will appear at the end of their respective chapters.) In 
this note, we shall present a set of curves giving spillover loss for several practical cases. 
For one feed case, we will also present a simple analytic expression which gives a good 
approximation to the spillover loss and is obtained by curve fitting. 

In spillover or directivity calculations, one must first determine the power radiated 

by the feed from the Poynting integral over a closed surface S, i.e., 

P = If s (ExH*)*ds (2- 1 ) 

If (E,H) were known exactly, the result of P would be unique no matter which surface S is 
used. In the Chu-Kirchhoff formulation [2, Eq. 11], the aperture field at the waveguide 
opening Si is approximated by the incident mode (e.g., TEn). The higher-order modes 
and reflected field are ignored. This is justified in the literature for circular waveguides 
with a>l X [2] and for rectangular waveguides [5]. Because of this approximation, the two 

power computations 

Pi = P calculated from (2.1) by using surface Si, and 

p„ s p eairnlareH from (2. 1) by using the infinite radiation surface S„ 

are generally different In fact 


1 Some of these results have been published in TREE Transactions on Antennas and Propagation, June, 
1990 [ 24 ]. 
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P x > P„ ( 2 - 2 > 

because Pi contains both P„ and the reflected power. In the limit afk -» Pi reduces to 
P„ as expected. In Figure 2.2, we plot PJPi as a function of afk for different waveguide 
feeds. In spillover calculations, we use P- instead of Pi as the input power so that the 
spillover loss goes to zero as 0 -> Jt. 3 

It should be noted that for field patterns the theoretical and experimental results 
diverge as the size of the aperture decreases. Agreement is good for waveguides as large as 
those for which results arc shown in Figures 2.3 to 2.7 [2], [5]. However, the results 
should not be extended to cover waveguides smaller than those for which results are 
shown. 

Spillover results arc presented in Figures 2.3 to 2.7 in the following manner. 


Table 2.1. Modes used in waveguides 


Fig. 

waveguide 

mode 

2.3 

circular 

TE U 

2.4 

circular 

TE 21 

2.5 

2:1 rectangular 

sum 

2.6 

square 

sum 

2.7 

4 square 

difference 


As an example, consider a circular waveguide for a monopulse feed, using the TEn mode 
and TE 2 i mode for sum and difference patterns. Let the reflector be of the symmetric 
parabolic type with f/D=0.4. The extended half-cone angle of die reflector is 64°. In order 


3 Reflector directivity calculations were sometimes carried out by using Pj (not P— ) as the feed input 
power. Strictly speaking this practice underestimates the directivity of the reflector antenna. For example, 
for a circular guide feed with a*0.7 X and excited by TE 21 , P 0 JP 1 read from Fig. 1 is 0.91, corresponding 
to a directivity underestimate of 0.4 dB. 




to keep the spillover loss less than 1 dB, die diameter of the circular guide should be at least 
0.90 X. This size waveguide feed would have less than 1 dB spillover loss for both the 

TEii and TE 21 modes. 

For the case of a circular waveguide excited by the TEn mode, an analytic 
expression has been derived by curve fitting. This expression is fairly accurate over the 
range where 20° < 0 < 70° and spillover loss is less than 4 dB. It reads 
spillover loss (in dB) = -a0~P 

where logioa = 3.78(a/X)“°- 629 (2.3) 

P = 2.50(a/A,)~°- 324 . 

Here 6 is the half-cone angle in degrees and a is the parameter defining the feed size. A 
comparison of the results obtained from the computer code and from the formula in (2.3) is 
shown in Figure 2.8 for several sizes of waveguides. 



Infinite sphere 
S 



Figure 2. 1 . The total input power radiated from the feed may be determined by 
integration over surface S„ or over surface Sj. 




Figure 2.2. Relationship b tween die power at infinity P~ and the power radiated through 
the end of the waveguide Pi for excitations of a circular feed and rectangular 
feeds. 







CIRCULAR FEED TE n 



Figure 2.3. Spillover loss for a circular waveguide excited by the TEn mode. 











SQUARE FEED TE 



Figure 2.6. Spillover loss for a square waveguide excited by the TEmmode. 
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(OP) uBAcmidS 


Figure 2.7. Spillover loss for a square difference waveguide. 
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Figure 2.8. Comparison of spillover calculated by numerical integration and that by the 
simple formula in (2.3) for the circular waveguide excited by the TEn mode. 
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3. REFLECTOR ANTENNA ANALYSIS 

In this chapter the method of solution for reflector antenna problems will be 
discuss^ After describing the problem and the various elements involved, the individual 
f- l f TTimts will be discussed one-by-one. This derivation is not original, having appeared in 
other sources [6], but it is useful for understanding the methods used in in our reflector 
antenna calculations. 

3.1. Description of Problem 

The geometry of the problem under consideration is shown in Figure 3.1. A 
reflector S is illuminated by the incident field from an array feed. The method of solution 
used is Aperture Integration (AI), which provides the same degree of accuracy as Physical 
Optics (PO) and avoids the caustics (infinite field in the main beam direction) that occur in 
the Geometric Theory of Diffraction (GTD). In PO, the induced surface currents on the 
reflector S are approximated by 

J,s2nxH*. C3- 1 ) 

These currents are then integrated to find the far field E‘(r)* In AI the reflected-field and 
diffracted-field contributions at a point P 2 on the planar aperture surface S a arc computed. 
This is done for points forming a grid over the aperture surface. A Fast Fourier Transform 
(FFT) is then used to obtain the far-field E*(r) [6-9]. AI has several advantages over PO. 
First of all, the use of an FFT (allowed by integration over the planar surface SJ is 
numerically efficient. The present formulation also allows for the use of multiple 
reflectors. In addition, AI can be used to obtain near-field information. Many of today's 
experimental measurements for reflector antennas are conducted in near-field ranges so this 
method can be used as an analytical check. Other advantages of the formulation of AI 
presented here are 

(i) The surface of the reflector may be arbitrary. 
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(ii) The edge of the reflector can be an arbitrary curve lying on an elliptical cone or 
cylinder. 

(iii) The divergence factor of the Geometric Optics (GO) field is correctly 
computed. This allows the feed to be placed away from the reflector dish 
where the divergence factor is not unity. 

(iv) The edge diffracted field is included here. Two uniform theories are used to 
keep the aperture field continuous from the lit to the shadow region. 

In Sections 3.2 to 3.6 the elements of the reflector problem will be examined. In 
the next section the source will be studied. This leads to the incident field H on the 
reflector. In Section 3.3 the method of describing the reflector(s) is put forth. This is 
necessary for obtaining the field at some point on the aperture surface. Finally an FFT is 
used to obtain the far field. 


3.2. The Reflector Surface 

The reflector surface is described by an analytical equation. This equation may take 
a variety of forms, depending on whether the surface is a type of conic or not In general, 
the surface is described by the equation 

z = do + d t x + day + d 3 xy + t^x 2 + d 5 y 2 + daiPU.y)] 47 ^ (3.2a) 


P(x,y) = s t + sax + sjy + s 4 xy + S5X 2 + s*y 2 + &jx 2 y + Sgxy 2 + S9X 3 + s 10 y 3 + 


s„x 2 y 2 + s 12 x 3 y + s 13 xy 3 + s 14 x 4 + s^y 4 


(3.2b) 


Some examples follow. 

3.2.1. Parabolic reflector 

In this case the reflector equation is given by 


z = C,+ 


x 2 + y 2 

4f L 


(3.3) 



where C a is a constant and ft is the focal length of the parabola. 
3.2.2. Spherical reflector 


16 


A spherical reflector may be described by 



where Q, is a constant and R is the radius of the sphere. 

3.23. Hyperbolic reflector 

This shape is commonly used for subreflectors. The reflector equation is 
z * wo + 

3.2.4. Conic section 

For a reflector cut from a conic section, the reflector surface is given by 




( 3 . 4 ) 


( 3 . 5 ) 


( 3 . 6 ) 


where e is the eccentricity and f is the focal length. 

3.2.5. Boundary 

Two types of boundaries are frequently used and receive special attention. In the 
first case, the boundary T is the intersection of surface S and an elliptical cylinder (Fig. 


3.2a). The parameters of the cylinder are 
(xc, yc) = center of the ellipse 
(Ki, K 2 ) = semiaxis along (x,y) direction. 

In the second case, the boundary f is the intersection of S and and elliptical cone (Fig. 


33b). The cone axis lies along the (y - y c ) - z plane and has the additional parameters 
(x = 0, y = 0, z = -p) = tip of cone 
63 ■ inclination angle of cone axis measured from z-axis 
(01, 62) = half-cone angles in the (x - Xc) - z and (y - y c ) - z planes. 



The boundary must be adequately defined in order to calculate the diffracted field ( see 
Section 3.5). 


3.3. The Source 

The source is assumed to have a well-defined phase center at point Pi that is the 
radiating point for a spherical wave (H 1 ^ 1 ). When an array feed is used, each radiating 
element must be considered separately. The scattered fields from each feed element are 
superimposed to obtain the total field at point P 2 in the aperture plane. It is desired to know 
the value of the wave (H*JE*) where it is incident on the reflector S. 

3.3.1. Feed source 

The surface current at the radiating aperture of the mth element of the feed array 
may be expressed as 

Jm(*,y) = Im(exac i¥ + e,b) ( 3 - 7 ) 

where (a,b,y) arc real and 

a 2 + b 2 = 1. < 3 - 8 ) 

The parameters (a,b,tp) are chosen to establish the feed polar iz a t ion. Table 3. 1 shows the 

values of (a,b,y) that correspond to commonly used polarizations. 


Ta ble 3.1. Various feed polarizatio ns 



a 

b 

Y 

linear x 

1 

0 

0 

lineary 

0 

1 

0 

RHCP 

1W2 

1/V2 

90P 

LHCP 

1W2 

\H2 

-90° 


The radiated electric field due to J m given by (3.7) is 
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Em * (3.9) 

where f m (0,<t)) is the active element pattern of the mth element The function f m may be 
approximately expressed by 


f m (0,<t>) = eeUeJ©) (ae^cos^ + bsin<|>) + e^tWe) (t*os<|> - ae jv sin<|)) (3 10 ) 

where 

Ueto( 0) = E-plane active pattern of the mth element 
UHm(0) = H-plane active pattern of the mth element 
There are two particular cases of Ueiu and Unm that arc of interest The first case is the q- 

feed [10], where die functions are approximated by 

UEmW) = (cosG)^- volt (3.1 la) 

IW0) = (cos0) qH - volt . (3. 1 lb) 

The second case is the numerical feed. Sometimes it is possible to measure the far-field 
pattern of an element In this case, if the pattern can be well approximated by its E-plane 
and H-plane patterns, it may be easier and/or more accurate to use measured data for UEm 
and Uan- In this case, the values for Unm and Unm are read from a table. The radiated far 
field of the feed array is the product of the element field (3.9) and the array factor, 

E«= X E m (r) I m ei k «u Pm 

m* 1 » v • / 

where Em is given by (3.9), I m is the complex excitation coefficient and 


e„ = sin0cos<t>e x + sin0sin<J>e y + cos0e 2 
Pm = X m ex + ym e j + Zme z . 


(3.13) 

(3.14) 


Another feed type that is commonly used is the open-ended waveguide. This feed 
is evaluated using PO. In this case, the incident fields in the waveguide are assumed to be 
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the dominant mode in the waveguide (e.g., TE1 1). The total field in the aperture is the 
incident field plus the reflected field. The far field is 
Er = 0 , 

.-i^r / i r* \/ 1 1 


E. = l + 1 (t^Xe) cos6 ] <n - cos * + N y sta * ) • 

E* = d ^t“ s9 + * (wX e ) ] (N « sin * ' N » cos * ) ' 


N is the vector 

N = (1 + D f Ej c ik<“k flc o*+y*“e«ii*) dS 


n h 6 
Ja 


(3.15a) 

(3.15b) 

(3.15c) 

(3.16) 


and t is defined by 


£L.„ ooe 

t for TE-modes and t = -s — 

coil Pmn 


for TM- modes. 


The phase constant is defined by 

Pmn = (k 2 -0 1/2 • (3.17) 

In this report the most-commonly used open-ended waveguide is a circular guide excited by 
a TE mode (Fig. 3.3). In this case the rectangular components of E 1 * are [7] 

E x = ^Y 2i [J m -i(K mn p) sin(m - lty + W^p) «“<*“ + . (3 lga) 

Ey= ^V .[J n> _ 1 ( 1Cmn p) COs(m- lty- InH-l^nmP) <»s(m+ ltyl lgb) 

where J m (x) is the mth-order Bessel function and TOwO = 0, where a is the radius of the 
waveguide. This leads die far-field expressions 


Eg = r' +X‘- t^ 056 )] ■ 

^ ^ /« 4 A 


(3.19a) 
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e» =r 1 [^ + COS0 - r(%i - cose)] 


Jm(*mn a > sinO) 
i _ f ^ sinO ^ 2 

'-y *mn ) 


cosm<|) . 


(3.19b) 


3.3.2. Power radiated 

We will assume for now a planar feed array. The total radiated time-averaged 
power of the array, assuming forward radiation only, is given by 


£ -2 

P rad = ^-J q 2 dojT E a (r) E^(r)r 2 sin0d<j> 


(3.20) 


with 



= 120tc ohms . 


Substituting (3.8) into (3.16) gives 


^rad 



(3.21) 


( 3 . 22 ) 


Defining the power radiated as follows, 
M M 

P«d = 

m=l n=l 


(3.23) 


a new term, the power matrix Amn, appears. 

3.3.3. Incident field on reflector 

Up until now, all positions have been determined in the feed coordinate system. 

However the reflector surface S is defined in the main coordinate system. The feed 

coordinate system can be related to the main coordinate system by an orthonormal 
* 

transformation matrix A 


[ej ej eJ] T = A[e x e, ej 7 


(3.24) 



Therefore a point P = (x,y,z) on reflector surface S can be expressed In the feed coordinate 
system by 
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[x F y F zV = Atfx-Xi) (y— yi) (z-z^f (3.25) 

where Pi ** (xi yi zi) is die location of the first feed element in the main coordinate system. 
Using (3.18), E* can be found at P. The H-field can then be found by the relation 

(3 26) 

The incident field is then converted from spherical coordinates to Cartesian coordinates. 
Finally, the H-field is converted to the main coordinate system by using 

=T 


[H x H y H z ] T = A [H^ Hy p 
This is repeated for each element in the feed array. 


(3.27) 


3.4. Geometrical Optics Field 

The GO field consists of two parts: the incident field H* and the reflected field H r . 
The incident field at the observation point P 2 on the aperture grid is taken to be zero, since 
it does not contribute to the far-field pattern. The first step for determining the reflected 
firlH is to locate the reflection point on the reflector surface S. 

3.4.1. Reflection point 

Given a starting point at the feed Pi = (xi yi zi) and an ending point P 2 on the 
aperture grid, a reflection point ■ (x y z=f(x,y)) on S may exist The vectors 

dj = e*(x - x x ) + e,(y - y^ + e 2 (z - z^ 

d 2 = e x (x 2 - x) + e,(y 2 - y) + ej(z 2 - f(x,y)) ( 3 .28) 

are the connecting vectors from Pi to Cf, and O r to P 2 , respectively. The condition on the 
reflection point is that the distance (di + 62 ) be a minimum, i.e.. 



(3.29) 


^•(d! + d 2 ) = 0 , J”(di + d 2 ) = 0 


which can be written explicitly as 

J-[ (x - Xl ) + [«x,y) - x,]|£} +^{(x-x 2 ) + [Kx,y) - xj]^} = 0 

(y - y i> + lf(x - y) - z > ] ^} + ^{ <J " y2 > + Wx ' y) " Zj, f'} = 0 • (3.30) 


The root of the two nonlinear equations in (3.30) gives the location of die reflection point 
For a given Pi and P 2 , there may be 0, 1, or more than 1 reflection point The system of 
equations in (3.30) could also be solved if the three points are collinear. An additional 
condition avoids this. 




Z - Zi z - z 2 
— + ~. 


>8 


(3.31) 


where 5 is a small positive number. 

A root may exist but be outside T. If 

l! £ 1 (3.32) 

then the root is inside the boundary T and is in fact a reflection point on the reflector. 
Otherwise the point is discarded. The parameter li is given by 




(3.33) 


if T lies on an elliptical cylinder and by 

li = I z | _1 {[(x - XckotOj ] 2 + [(y - y c )cot 0 2 ] 2 }T ( 3 . 34 ) 

if T lies on an elliptical cone. The parameters Xc, y c , Ki, K 2 , 0i, and 0 2 describe the 


boundary. 

3.4.2. Formula for reflected field 

The reflected magnetic field at P 2 is given by 
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H r (P 2 ) = (DF)e _jlcdj (H*(O r ) - 2[H‘(O r )-eN]e N } (3.35) 

where DF is the divergence factor and eN is the surface normal of the reflector at Cr. The 
normal is chosen to be pointing towards the source so that (eN«i) is always greater than 
zero. Explicitly, es is given by 

eN = A(-f x e x - f y e y + z) (3.36) 

where 

1 

A = + (fjf + fJ+l) 2 (3.37) 

and the subscript x of f x , for example, means partial derivative with respect to x. The 
divergence factor in (3.35) describes the spreading out of the wavefront and is expressed as 

™ I 1 

l + (d 2 /R r 2 ) (3.38) 

where the square roots take positive real, negative imaginary, or zero values (so that DF is 
positive real, positive imaginary, or infinite). (R r i, R r z) the principal radii of curvature 
of the reflected wavefront passing through O r . Their computation is covered in Section 

3.4.3. 

3.4.3. Curvatures of reflected wavefront 

The formulas from [1 1] are used in calculating (R r i, R r 2). The three orthonormal 
base vectors of the incident pencil are chosen to be (Figure 3.4) 

I X Cj(z - 2[)“ti(X”Xi) 

** |e y xej| [(z-z 1 ) 2 + (x-x 1 ) 2 ) 1 ' 2 (3.39a) 

, _ e^xe} 

I el x ej I (3.39b) 

, _ e^x - xi) + e T (y - y t ) + e,(z - z t ) 

3 [(x-x 1 ) 2 + (y-yi) 2 + (z-z 1 ) 2 ] 1/2 



(3.39c) 
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where (x, y, z) are the coordinates of the reflection point 0 T . Those of the reflected pencil 
are chosen to be 


ej = - 2(e| e N )e N 

(3.40a) 

e 2 * *5 - 2(©l e N )e N 

(3.40b) 

r *x( x 2 “ ») + e y^2 " y) + ^ 

I(x 2 -x) 2 + (y 2 - y) 2 + (*2 " z ) 2 ] 1/2 

(3.40c) 


Note that (3.40) establishes a left-hand system, i.e., 

*1 * *2 = ”®3 * 

This does not affect the final solution of (R r i, R r 2 >- The three orthonormal basis 
vectors of reflector S at (f are chosen to be 

S _ 

ei ”(l + ^) 1/2 (3.41a) 



s s 

*3 x e l 



e| = e N 

From (3.39) and (3.41), the elements 

pL = • mji = 1, 2, 3 

can be calculated with the results 

. (z - z t ) - f x (x - Xx) 

Pll “ 0+^) ul I(*-x,) 1 +(*-*i) S ] 1 ' 2 

j -fy[(x - X t ) + f x (z - Z t )] 

Pl2 ” A(1 + $ 1/2 [(x - Xl ) 2 + (z - z,) 2 ] 1/2 


(3.41b) 

(3.41c) 


(3.42) 


(3.43a) 


(3.43b) 
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” (x - x x )(y - yi) - f x (y ~ yiKz ~ z i> 


P 21 - 


(1 + # l/2 {(x - x x ) 2 (y - y x ) 2 + Kz - zi) 2 + (*- x x ) 2 ] 2 + (y - yi> < z “ z i> > 


^ 1/2 


(3.43c) 


P 22 = 


f,f v (x - x,)(y - y,) +0 + f,)[fa - n) 2 + (» - ill) 2 ! - fyfr - ?!><* - *1> 

A(1 + # 10 {(x - x,) 2 (y - y,) 2 + Kz - z,) 2 + U - *,)¥ + <y - yi»*<* ' 

(3.43d) 

P 33 = — A[f x (x - x x ) + - y x ) - (z - z x )] • (3.43e) 

The first four elements ((3.43a ) - (3.43d)) form the 2 x 2 matrix P\ Because of the 
particular choice in (3.40), we have P = P. The curvature matrix of the incident pencil is 


Q = 


dl 1 0 

j-i 


L 0 d XJ 

and the curvature matrix of reflector S at C? is 

_S cA 2 (G - F 2 /E) A(eF - fE)/E 
Q ” A(eF - fE)/E A 2 (gE - 2fF + ef^/E) 


(3.44) 




where 

E = l + £, F = f x f y , G-l+f 2 

eS-Afjm , f=-Af xy , g = Afyy 

The desired curvature matrix Q 1 may be calculated from the following matrix equation 

Q = Q + 2 p 3 3 [( P )V Q ( P r 1 • (3.46) 

The four elements of $ may be denoted by 
[Q11 Q12I 

Q LQ21 Q22J ( 3 . 47 ) 

Then the desired radii of curvature of the reflected wavefront at Cf arc given by 

— , — = 4 { (Q11 + Q22) ± V ^Qn + 022^ 2 ” ^QllQ22 “ Q12Q21)} * 

R r x R 2 2 (3.48) 
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Both R r i and R r 2 arc real. If the radius is positive, then the corresponding normal section 
of the reflected wavefront is divergent. If the radius is negative, then the corresponding 
normal section of the reflected wavefront is convergent 

If there is more than one reflection point, then the total reflected field is the 
superposition of the reflected fields for each reflection point If no reflection point exists, 
then the reflected field is zero. If the reflected point is near the boundary, then the reflected 
field is calcu lated as normal, but the diffracted field is adjusted later by using UAT so that 
tiie total field is correct 


3.5. The Edge-Diffracted Field 

In addition to a reflected field at P 2 , there is also an edge-diffracted field. The first 
step for finding the diffracted field is to find the diffraction points on the boundary T. 


3.5.1. Diffraction points 

Consider a source point Pi at (x'i, y'i, z'i) and an observation point P 2 at 
(x' 2 , y’ 2 , z' 2 ) with coordinates given in the prime coordinate system (Fig. 3.5). A 
diffraction point O d with coordinates (x\ y\ z') can be determined from Fermat’s 
principle, much as the reflection point was 




(3.49) 


where d 3 and dLj are distances from Pi to O and 0 to P 2 , respectively. This may be 
rewritten as 


81-*! 

d 3 


+ 


gr-xiVgi f 82 - yi 

d 4 J<5<r l d 3 


g2-y2l d g2 . ( g3“ z l 

d 4 JW l d 3 


+ 


g3” *2^3 

d 4 JW 


= 0 
(3.50) 


where (gj, g 2 , g 3 ) are parametric functions describing x\ y\ and z' in terms of 4>'. A root 
of this equation is a diffraction point Depending on the geometry, as many as four 
diffraction points can exist 
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3.5.2. Formula for diffracted field 

Each diffraction point has associated with it a contribution to die diffracted field n. 


The formula is 




(3.51) 


where g is a cylindrical wave factor 

+ • (3.52) 

The other factors are explained in more detail in the following sections. 

3.5.3. Divergence factor 

The square root in (3.51) can be positive real, negative imaginary or zero. Ri is a 
radius of curvature of the diffracted wavefront passing through O d and is found from the 
equation [11] 


_1 1 _ 

Rj d 3 


sin p 


(3.53) 


where P is the angle between tangent t and d* The term k is the curvature of curve T and 
e B is the normal of this curve at O d . This leads to a final expression for the divergence 
factor of the diffracted field 


DFD = 


yTTwvRT) 


(3.54) 
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where 


c= £i+ * J [g 3 (g,s; - g.gi> - g 2 <gig2 - g2g".)i 

d 3 P*(«tap) 2 ^ d 3 d * > 


+ ^ ■ y * + — J“jlgl(glg2 “ tllO * gj<g2g3 - g3*2>l 
+ ^j^jlgi(gigj - g3g2> ~ 8i ( 83gi _ *1*3» 
where a g-prime (g 1 ) represents a derivative of g with respect to $ , 

n i 


(3.55) 


p = V tei ) 2 + te 2^ 2 + ( fo ) 2 • 


(3.56) 


(3.57) 


3.5.4. Diffraction coefficients 

The hard and soft Keller's diffraction coefficients D h and D* are defined in Eq. 


( 5 . 22 ) of [ 11 ] as 


D ,,h = x* + X* = “^cc^te ” sec^te + $ ) • (3.58) 

The angles <j>and<|> l are shown in Figure 3.6. They are calculated from the relations 


(3.58) 


cos 4 >' = (- 1 ) 


(Proj d 3 )-(e t x e n ) 
| Proj d 3 1 


(3.59) 


(Proj d)-(e t x e B ) 

cos * ikjdi 


(3.60) 


Here Proj d 3 is the projection of d 3 on the plane perpendicular to e». This can be rewritten 


2 , c2 . o2 1 


COS^ 55 [S 31 M 1 + S 32 M 2 + S 33 M 3 ][S 3 i + S 32 + S 33 ] 


(3.61) 
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COS()> = [S41M1 + S42M2 + S43M3HS41 + S42 + S43] 


(3.62) 



M 1 = |rt N 2g3 ~ N 3g2^ 

M 2 = ^[N3g' 1 -N lg ;] 

M 3 = |-[N 1 g 2 -N 2 g’ 1 ] 

and Ni, N 2 , and N 3 are the components of eN, the normal to surface S, in the (x\ y\ z') 
coordinate system. 

The solutions of <|> and (p 1 are subject to the following tests: 

(i) 0 < <|> l < Jt if T 3 £ 0, and it < tp 1 < 2ti if T 3 < 0, where 

T 3 = (-Proj d 3 )-e N = -S 31 N! - S 32 N 2 - S 33 N 3 . ( 3 .63) 


(ii) 0 < <(> < Jt if T4 ^ 0, and 7 t < (j) < 2 ti if T4 < 0, where 

T 4 = (-Proj d 4 )-e N = S 41 N { + S 42 N 2 + S 43 N 3 . ( 3 .64) 

(iii) If the observation point P 2 is exactly on the incident shadow boundary, then 

<p - <p* - 7C = 0 . (3.65) 

(iv) If the observation point P 2 is exactly on the reflected shadow boundary, then 

$ + f - jc = 0 . (3-66) 
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When P2 is on the incident shadow boundary (case iii), the factor x‘ becomes 
infinite. The resulting computational problem is avoided by moving P2 slightly whenever 
this occurs. Similarly, % r becomes infinite in case iv, and the same procedure eliminates 

die problem here also. UAT wiU be used later to modify the diffracted field H* to H . 

3.5.5. Spherical components of incident field 

Fields H l p and H‘ a from (3.51) are the two spherical components of the incident 

magnetic field H 1 evaluated at O d in the directions fi 1 and a 1 . The base vectors are shown 
in Figure 3.7 and may be calculated from 

e p ‘= (e 0 'xdj)/d 3 (3.67a) 

e a » = (e t x d 3 )/(d 3 sinp) (3.67b) 

Then it can be shown that 

H P = + H2g2 + H3g3] (3.68a) 


H “ = Pd 3 W {Hl[g2(Z " Zl) ' g3(y " yi)1 + H2lg3(x “ Xl) ~ gl(z " Zl>1 

+ H 3 [gi(y - y’i) - gi(x’ - *i)] 


(3.68b) 


where (Hi, H2, H 3 ) are the components of H 1 in the primed coordinate system. 

3.5.6. Rectangular components of diffracted field 

The only part that needs to be determined now are the spherical base vectors for the 
diffracted field, ep and e<x. They are defined by equations similar to (3.67), namely, 

ep * (e a ‘ x d4)/d 4 (3.69a) 

e a = (e t x d 4 )/(d 4 sinP) . (3.69b) 

Once these vectors are determined, they must be expressed in terms of the rectangular 
components (e t , e, f e*) so that the diffracted filed may be superimposed upon the 
reflected field. The corresponding diffracted E-field is then 


Ed(p 2 ) = ^[ Hd(p 4 ) x d 4 ] . 


(3.70) 
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3.5.7. Detour parameter 

As mentioned in Section 3.5.4, when the observation point is near the reflected 
shadow boundary, the diffracted field H d is not valid. The "detour parameter" is used to 

determine if this happens. From Section VI of [1 1], 

% = e[k(d 3 + <14 - d! - d 2 )] 1/2 . (3.7 1) 

Here e is the shadow indicator of the reflected field and is +1 or -1 if the point is in the 


shadow or lit region of the reflected field, respectively. It may be shown that 

e = Sgn[-cos±(* + 1i i >] (3 . 72) 

When the caustic of the reflected field falls on the reflected field shadow boundary, ^ is 
imaginary; otherwise, it is real. Following die numerical study in [12], 

Ul*2 (3-73) 

is the dividing line. This means that if I ^ I >2, then the field H d is valid. If not then H 
must be replaced by H^. When ^ is small, (di + d2) and (d 3 + <14) from (3.71) are nearly 
equal and an alternate formula for £ is used, namely. 


1 ; . V ^[1 + ( ( V R i>l t 

% - -cos-r-(<{> + $ ) sin(S V2k - — . =r 4 

V 1 + W2/R1) V 1 + (d2/R T 2> 


0 , 


(3.74) 


where the square root in the numerator is the divergence factor of the diffracted field and 
the square root in the denominator is die divergence factor of the reflected field. 

3.5.8. Uniform asymptotic theory 

We shall calculate by the UAT developed in [13-15], namely, 

H d (P 2 ) = H d (P 2 ) + [f(£) - H£) - ^(1 - e)]H r (P 2 ) , 


where F is Fresnel integral defined by 

F(z) = 7t~ 1/2 e^ 4 f* c~ }t dt and 
Jz 


(3.76) 



32 


• (3.77) 

The factor (1 - e)/2 in (3.75) is one if P 2 is in the lit region and zero if it is in the shadow. 
The multiplier for H r goes to zero as £ goes to infinity, so that H° becomes equal to H d . 
Near the boundary, H D is finite and compensates to make die total field continuous. 


3.6. Secondary Pattern Computation 

From field equivalence principles, solutions for the far field may be obtained if the 
tangential fields E a and H a on the aperture plane are known. This can be done directly. 
The following vector quantities may be defined 

f(«,v>= ff E.(x.y) e ik(l “ +vy) dxdy 

JJ s. (3.78) 

*(u,v)- ff H t (x,y) c* k(ra+v>1 dxdy 

JJ S. (3.79) 


where 

u = sin0 cos<(> 
v = sin0 sin<|> 
k = 2 ji/X 

0 , <p = sperical coordinates of the far-field point 
Since the aperture fields are tangential and the far field is assumed to be in the z-direction, f 
and g will have only x- and y-components. There are three equivalence principles and each 
iwrft to a different pair of equations for the the far-field E-field components Ee and E*. 


(1) Using E a and H. 

-jkr 

Ee = J 4 ■ [f x cos<p + fySin<t> + Z^os© (g y cos4> - g x sin<t»)] 

-jkr 

J 4jn . [cos0(fycos<|> - f x sin<f>) - Z<, (gySin<f> + g x cos<|>)] 


(3.80a) 

(3.80b) 


(2) Using H, 


E 

E ®” 4jtr 


Z^cosO (gyCos(f> - g x sin<>) 


(3.81a) 



33 


ike - ^ 1 

E*= J 4jn . Z;(g y sin<j> + g x cos<l>) 

(3.81b) 

(3) Using E„ 


jke _ikr 


E e = (f x cos(|> + f y sin<J>) 

(3.82a) 

]kc -jkr 

E^ = cos0(fyCos<|> - f x sin<t>) . 

(3.82b) 


This method is exact if the aperture fields are known everywhere. In practice the aperture 
plane must be truncated in order to employ the FFT. 

3.6.1 Fast Fourier transform 


Some manip ulation of the integrals in (3.78) and (3.79) must be done in order to 
use the FFT. When the components of the aperture field have been correctly scaled, the 
fif ths ran be approximated by a double summation over die grid points in the aperture 
pin nr. The grid is two-dimensional and covers a portion of the aperture plane that is 
slightly larger than the reflector, i.e., the grid extends past the projection of the boundary T 
on the aperture plane. Each term includes a scalar constant factor, the Fourier coefficient 
Qnn. After some more manipuladon and evaluation of the integrals, the expression for f* is 


M 2 

f„(u,v) - K X E C-'*”"’ 

i»bN| ro=Mj 


sin y(mX+u(x 2 -x 1 ))jsin^-^(nX+v(y 2 -yi)) 
7i(mX.+u(x 2 -x 1 )) Jt(nX+v(y 2 -yi)) 

X X 

(3.83) 


where 

K = (x 2 - x)(y 2 - *> . (3 . 84 ) 

S imilar results can be obtained for the other far-field components. 

3.6.2. Polarization of secondary pattern 

The polarization of the secondary pattern is described by the elements (a p ,b p ,t|/ p ). 
hi the planar aperture of the main reflector at z = 0, the tangential electrical field polarization 


can be written 
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E ap (z=0) = (e 1 a p e jvP + e,b p ) g(x,y) , ' (3.85) 

where g(x,y) is the aperture field distribution in the xy plane. From [16] we can define the 
orthonormal reference and cross-polarization vectors of the secondary pattern by 


eji = e^aV^cosf + l/sinf) + e^-aV^sin^’ + bPcosf) (3.86) 

«c = ^(aV^rinsf - t^cos*’) + e<,(a p e -jyP cos<t>' + t^sinf ) . (3.87) 

The polarization parameters (a p ,b p ,y p ) of the secondary pattern are related to the 
parameters (a,b,y) of the feed by the relationship 

a p = a, b 55 = b, y p = \|/-Mi • (3.88) 

For example, this means that an RHCP feed produces an LHCP secondary pattern. This 
relationship is true only when a single reflector is involved. For the multirefiector case this 
may change. It is important to keep track of the number of reflections in the system in 
order for the secondary pattern polarization to be correct 


3.6.3. Directivity 

Once the ref-pol of the secondary pattern for a single element has been found, the 
process can be repeated for all of the elements. The total ref-pol of the secondary pattern is 


M 


Er = ^ Enn,!,,, 


ro=l 


(3.89) 


where Erbi is die ref-pol field for the mth element assuming excitation In, = 1, and In, is the 


actual feed excitation. 

The directivity of the reference polarization is defined by 
_ 4 ic(Er) 2 

D^e,*)- 2^ (3 .90) 

where P r ad is the power radiated by the feed. It can be calculated by a brute-force 
integration over either half-space or all space. The cross-pol directivity may be calculated 
by a similar method. 
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3.7. Concluding Remarks 

The main usefulness of the AI method is that it is very general in nature. Until 
recently, when the currents on a reflector were found, steps would be taken to reduce the 
resulting integral so that its computation would be possible. These steps were usually 
b ased upon some aspect of the reflector geometry. With the advent of the supercomputer, 
this is no longer necessary. A brute-force FFT can be used to evaluate the fields in the 
aperture without any geometry-based simplifications being made. Since the FFT is not 
d epe ndent upon the reflector shape, this method may be applied to arbitrary reflectors. In 
addition, the method may be repeated in order to include additional reflectors. 
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Feed Array 


Figure 3.1. Secondary pattern of a reflector antenna using Aperture Integration. 
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Figure 3.3. Coordinate system used for open-ended circular waveguide. 




41 



Figure 3.6. Projection of Figure 3.5 on a plane perpendicular to tangent e t . 
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Figure 3.7. Spherical base vectors. 
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4. WIDE-ANGLE SCANNING FOR REFLECTOR ANTENNAS 

Traditionally, reflector antennas are designed for limited scan 4 . A symmetrical 
parabolic reflector with f/D=*0.4 can only scan ±5 beamwidths (BW) with less than 2 dB 
loss [4]. If the reflector diameter is 1,000 X, the ±5 BW scan corresponds to only ±0.5°, 

which is a very narrow field of view. 

In some future applications, the antenna requirements will be quite different from 
what they presently are. One example is the NASA Earth Science Geostationary Platform 
Project The preliminary antenna specifications are as follows: 

Frequency range 20 GHz - 200 GHz 

Antenna diameter 15 m (1,000 X - 10,000 X) 

Scan range 1 i2° (133 BW - 1330 BW) 

Scan range 2 18° (1133 BW - 11,333 BW) 

Note that the scan requirement has been significantly increased from the traditional value of 
15 BW. Usually a phased array design is used to satisfy specifications such as these. 

A phased array antenna design is an order of magnitude more complicated than a 
reflector design. This is due to the large number of array elements and die beam-forming 
network con tained in die design. Reflector antennas have the additional advantage of being 
less expensive and lighter in weight than phased arrays. Therefore, it is desirable to use a 
reflector antenna design if at all possible. The question then is "Can a reflector antenna be 
desig ned that is capable of meeting these specifications?" 

This paper examines and compares six different reflector designs. We intend to 
show how far die reflector performance can be stretched. The object is to achieve a wide- 
angle scan that will satisfy requirements such as those listed above. The first three designs, 
PI, P2, and P3, are parabolic single reflector designs. These three designs are considered 

4 Some of these results have been published in Microwave and Optical Technology Letters in July 1990 
[36], 
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in Section 4.1. The first design is a center-fed, single-element feed design with f/D=2 
(Fig. 4.1a). Scanning is accomplished by mechanically tilting the reflector. The second 
design has f/D=l and uses a 19-element cluster feed but, otherwise, is similar to the first 
design (Fig. 4. lb). The third design is an off-set reflector with flD=2 and an electronically 
scanned cluster feed (Fig. 4.1c). 

The last three designs. Cl, C2, and C3, are dual reflector cassegrainian designs. 
Urey are considered in Section 3.2. The three designs all use the same reflector geometry. 
The first design scans by mechanically tilting the main reflector (Fig. 4.2a). The second 
design scans by mechanically tilting the subreflector (Fig. 4.2b). The last design scans by 

tilting both the main reflector and the subreflector (Fig. 4.2c). 

Data on extremely wide-angle scans of reflector antennas are scarce in the literature. 
Hung and Mittra [17] in 1986 did analyze a center-fed symmetrical parabolic reflector with 
a cluster feed and calculate patterns up to a hundred-beamwidth scan. We have verified our 
single reflector computer code by comparing our results with theirs. 

4.1. Single Reflector Antennas 

4.1.1. PI: symmetric parabolic reflector with flD=2 

Two contributing factors to poor scanning ability are (i) short focal length and (ii) 
high offset For these reasons, the first design considered is a symmetrical parabolic 
reflector with an unusually long focal length. Design PI has a diameter D= 1,000 X and a 
focal length f=2,000 X, thus giving f/D=2. The feed is a long circular open-ended 
waveguide with radius a=3 X. A study was done of the directivities and beam efficiencies 
corresponding to various feed radii The results are shown in Figure 4.3, with directivities 
converted to antenna efficiency. Antenna efficiency is defined as the fraction of the 
nominal directivity that the given directivity is, namely, 
r\ mt as (Directivity )/(tcD/X) 2 


(4.1) 
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In this case, the nominal directivity (xD fk)^ is 69.9 dB. The radius value chosen was that 
which maximized beam efficiency. The antenna has a half-power be am width 
HPBW=0.07°. Beam efficiency is calculated as the fraction of power hitting the reflector 
that is contained in the beam defined as being 2.5 times as large as die HPBW [18]. In this 
case, the beam has a half angle of approximately 0.09°. Note that this definition of beam 
efficiency does not take into account spillover loss. The chosen radius value of 3A, 
produces the highest beam efficiency, T|=0.91. Scanning is accomplished by tilting the 
main reflector. Tilting the main reflector by a degrees results in a scan angle 0*2a. The 
main advantage of tilting the reflector instead of moving the feed element is that the scan 
angle is twice die angle of the tilt If the feed were moved, then the angle of scan would be 
equal to the angle that die feed was moved through. This is referred to as the mirror effect 
Since for any reflector design, the scan loss increases as the feed moves away from the 
reflector's focal point a significant reduction in scan loss is gained by tilting the reflector 
instead of shifting the feed. 

Features of reflector PI are 

1 . Virtually no feed blockage due to the small size of the single element 
feed. 

2. Depending on the exact arrangement there is a lossy transmission 
distance between the feed and the receiver/transmitter. To avoid 
excessive transmission loss at high frequency applications (60 GHz or 
more), it may be necessary to connect the feed and the 
receiver/transmitter via a beam waveguide. 

3. Because of the mirror effect the scan range is twice as far as for the 
conventional shifted feed design. 

The radiation pattern for the on-axis beam is shown in Figure 4.5. The radiation 
pattern is calculated by a standard physical optics reflector code [19]. The directivity is 
66.7 dB which includes the following losses: 
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Nominal directivity (xD/X) 2 

69.9 dB 

Feed spillover loss 

- 0.6 dB 

Amplitude taper over reflector surface 

- 2.6 dB 

Directivity 

66.7 dB 


The above directivity, as usual, does not include the loss due to the feed transmission line. 

The 3 X radius feed produces a pattern that has a null before the edge of the reflector (Fig. 

4.4). This pattern results in a sidelobe level of -31 dB. It is a commonly used rule-of- 
thumb that to maximize beam efficiency, the first feed pattern null should lie on the reflector 
edge, which is at 0 max =14.25°. The first null lay on the reflector edge for feed radius 

a-2.3 X. This value produces close to a maximum in beam efficiency (see Fig. 4.3). 

This reflector has extremely good scan characteristics because of the long f/D and 
the mirror effect The scan loss is only 0.6 dB at 00=8° (Fig. 4.6), corresponding to a 1 14 

beam width scan. The sidelobe level does increase from -31 dB to -13 dB as expected. At 
a larger scan angle 0 o =2O° (286 beam width scan), the scan loss is 5.1 dB and the pattern is 

badly distorted with a high shoulder (Fig. 4.7). 

4.1.2. P2: symmetric parabolic reflector with f/D=l 

A drawback of PI is its excessively long focal length (2,000 X). Now let us reduce 
it by one half, giving a ffD*l. Then wide-angle scan is possible only if a cluster feed is 
used. A brief explanation of the cluster feed concept is in order at this point The feed 
cluster consists of N identical elements with complex excitations 

I*[Ii, I 2 I N 1 (4 * 2) 

We wish to determine I so that, when the beam position is at 0 O , a prescribed 
antenna parameter such as directivity, beam efficiency, or sidelobe level is optimized. To 
this end, let us introduce an element secondary pattern vector E such that 


E(0 o HE 1 (0 o ), E 2 (0 o ),...,E N (0o)] 


(4.3) 



where E2 (0q), for example, is the co-polarization secondary pattern in direction 0 O when 
element 2 is excited with 


I 2 =l 

^ = 0 , for all m * 2 . 


There exist three methods for determining I in the literature. 

Cl) rnniugate Field Matching [20-27]. The cluster excitation is simply set equal 
to the complex conjugate of E(0 q), i.e., 

I=[E(0 o )]* (45) 

Strictly speaking, such a choice of cluster excitation does not optimize any 
particular antenna parameter. For practical purposes, however, it does lead to 
nearly optimum directivity in most cases. 

Cri) Optimum Directivity [28]. For a feed cluster with prescribed primary 
patterns and element locations, the directivity in direction 0 O is optimized by 

choosing 

1 = A' 1 [E (0 O )]* (4-6) 

where A is a N x N square matrix with elements 

(4-7) 

where C is a normalization constant, and the integration is over 4it - radiation 
sphere. When the element spacing of the cluster is large (a few wavelengths), 
matrix A is nearly an identity matrix. Then the solution in (4.6) reduces to 
that in (4.5). 


(iii) SiHrlnhe Control [29,30]. The element secondary pattern vector E in (4.3) is 
normally calculated in a transmitting approach. By reciprocity, it can be 
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calculated equally as well in a receiving approach when the reflector is 
illuminated by an incident plane wave from direction 0 O . In the receiving 

approach, there is an additional advantage that the amplitude of the plane wave 
can be tapered. It is found that the amount of taper controls the sidelobe level 
of the final secondary pattern when the whole cluster is turned on. 

Return now to P2 in Figure 4.1. A 19-element cluster feed is used. The individual 
elements are circular feeds with radius a=l .2 X. This value is chosen to maximize 
directivity for a single feed scanned on-axis. This radius feed also produces a relatively 
good beam efficiency with q=0.88. The maximum beam efficiency was 11=0.89, which is 
recorded for a feed with radius a=1.3X. Sidelobes for the a=1.2 X feed were -23 dB. This 
was not nearly the best possible sidelobe value, as a feed with radius a=1.5 X had sidelobes 
that were -32 dB, along with q=0.87. The primary pattern of the a=1.2 X circular 
waveguide feed was approximated by a (cos0)9 pattern with q=9.5. This value of q gave 
good sidelobe matching but the main lobe was 0.6 dB higher, with a maximum directivity 
of 68.9 dB. Spillover loss for the cos^ feed pattern was 0.4 dB for a=1.2 X. The cluster 
feed is used to help compensate for the higher scan losses that result from the lower f/D. 
At small scan angles only the center feed element has a relatively strong excitation (Fig. 
4.8a). For an 8° scan (i.e., the reflector is tilted 4°) only two of the outer ring elements 
have sig nifi cant excitations (Fig. 4.8b). This indicates that for scans under 8 , a 7 -element 
cluster feed would probably work almost as well as a 19-element feed. When the reflector 
is tilted 10° for a scan angle of 20°, nearly all of the elements are excited (Fig. 4.8c). At 
6 0 =8°, the scan loss is 3.7 dB (Figs. 4.9, 4.10), and at 0 o =2O°, the scan loss is 7.4 dB 

(Fig. 4.11). Very similar scan loss results were obtained with a=1.5 X feed. This is not as 
good as the results for PI, but it is only a few decibels worse. The advantage of P2 over 
PI is that the focal length has been cut in half. The disadvantage is that a 19-element feed 
is much more complicated than a single element feed. For both of these center-fed designs 
the feed blockage is negligible. 



Though design P2 has a higher scan loss at 0 o =2O° than PI, the beam is less 
distorted (see Fig. 4.7 and Fig. 4. 1 1), because at scan angles of this size, the cluster feed is 
able to form a much better beam pattern than does a single element feed. For angles below 
8 0 =8°, there is no benefit to design PI from using a cluster feed. However, scan loss 

could be reduced for large scan angles by using a cluster feed. 

4.1.3. P3: offset parabolic reflector 

Design P3 is an offset parabolic reflector. Offset height must be kept as small as 
possible to avoid intolerably high scan loss. Unfortunately, small offset leads to serious 
feed blockage. A possible way out of this dilemma is to use two identical reflector 
antennas: one for scanning up and one for scanning down as sketched in Fig. 4.1. The 
focal length is 2,000 X and the reflector diameter is 1,000 X, for a ffd= 2. The offset height 
is zero. In contrast to PI and P2, this design utilizes electronic scanning. This means that 
a large feed array is used. Up to 19 elements are excited at any time. In order to cover a 
scan range from 0 o =O° to 9 0 =8°, a semicircular array with a radius of 283 X must be used. 

The individual feed element has a radius a=1.065 X, meaning that roughly 65,000 elements 
are needed for the entire device. This feed size is chosen so that if the excited element is 
turned off and an adjacent element is turned on, then the beam is scanned 1 BW. This 
antenna has an on-axis directivity of 66.2 dB when a single element is turned on. The 
spillover loss is 3.5 dB. This is quite high since the feed element is so small. Note that 
this antenna has roughly the same f/D ratio as PI, which uses a feed that is three times 
larger in radius. 

The advantage of the electronic scan is well-known: it is fast and inertialess. 
However, electronic scanning forces the use of a feed array that has half the diameter of 
one of the reflectors used. Therefore, this design uses about a fourth as many elements as 
a phased array with the same aperture size. The savings in complexity are almost lost. In 
addition, the overall volume occupied by this antenna is much larger than that needed by the 
previous designs. For the on-axis beam, only the center element of the 19-element cluster 
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is significantly excited (Figs. 4.12, 4.13), with a directivity of 67.3 dB. This is slightly 
higher than that excited by a single clement feed (66.2 dB). Although excited with small 
excitations, the surrounding elements do help to reduce the spillover (Fig. 4.15). 

Another problem is that the design puts a physical limitation on the maximum 
scanning angle, because the feed elements do not move. The previous designs could have 
been scanned farther than 20° if it had been desired. This design has a directivity of 67.3 
dB, with BW=0.06° and sidelobes at -18 dB. Patterns were computed at scan angles 
0 o =O°and 0o=8°. The feed excitations used to obtain these results are shown in Figure 
4.12. At 00=8° the scan loss is already 6.3 dB (Fig. 4.13, 4.14). The advantages of 
electronic scanning are that it is quicker than mechanical scanning and that it will not upset 
the equilibrium of the spacecraft since there is no physical motion. Some of the drawbacks 
listed above could be avoided by mechanically moving a 19-element feed cluster instead of 
electronically scanning. However, this design has much more scan loss at 00=8° than 

designs PI and P2. 

4.2 Dual Reflector Antennas 

4.2.1. Cl: cassegrainian reflector with tilted main reflector 

Design Cl is a dual-reflector cassegrainian antenna. The main reflector is parabolic 
with a focal length of 2,000 X and a diameter of 1,000 X for a f/D=2. A Cassegrain 
antenna may be considered as a folded version of a parabolic reflector. In many 
applications, it is desirable to reduce the length of the antenna and to place the feed directly 
behind the vertex of the main reflector. These are the reasons for folding the antenna. 
With f/D=2 for the present case, it is not possible to fold the feed close to the vertex 
without either excessive spillover loss or an excessively large subreflector or even both. In 
the present design (Fig. 4.2), the hyperbolic subreflector has a diameter of 1 15 X and is 
located 1,650 X from the main reflector vertex. The circular feed has a radius a=1.5 X and 
is located 1,300 X from the main reflector vertex. This feed size is chosen to produce 10 



dB edge taper on the subreflector. Directivity for this design is 67.1 dB, with BW=0.06° 

and a -18 dB sidelobe level. Scanning is accomplished by tilting the main reflector. The 
scan angle Qg-aj, where (Xj is the angle at which the main reflector is tilted. The 

performance of Cl, shown in Figures 4.15 and 4.16, is similar to that of PI, the unfolded 
version of Cl. The use of the subreflector does change the aperture taper. Consequently, 
the sidelobes of Cl and PI are different 

4.2.2. C2: cassegrainian reflector with tilted subreflector 

Design C2 has the same geometry as Cl. Scanning is accomplished by tilting the 
subreflector instead of the main reflector. Due to the substantial difference in size, tilting 
the subreflector is much easier mechanically than tilting the main reflector. Electrically, 
however, tilting the subreflector for wide-angle scan is not feasible because 

(i) the subreflector must be tilted by a much larger angle otj in order for the beam 

to scan. The approximate relation between the two angles is 

0O-OC2/M 

where M=D main /D sub =magnification factor. 

(ii) When the subreflector is tilted by a large angle, there is an excessive spillover 
loss. 

In this case M=8.7. The scan loss is quite high. At 0 O =1.75°, the scan loss is 6.6 dB 
(Fig. 4.18). At 00=3.32°, the scan loss is 36.3 dB. This would seem to indicate that 

tilting the subreflector is not a viable option for wide-angle scan. 

4.2.3. C3: cassegrainian reflector with both reflectors tilted 

Design C3 has the same geometry as Cl and C2. Scanning is accomplished by 
tilting both the subreflector and the main reflector. The idea is to use the main reflector for 

coarse scanning and the subreflector tilting for local scanning within a small angular region. 
The scan angle 0 o *2(aj+O2/M). Given aj and o^, the scan loss can be obtained by 

looking at the results for Cl and C2. 
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4.3. Conclusions 

We have studied the wide-angle scan ability of the six reflector antennas shown in 
Figures 4.1 and 4.2. All reflectors have a circular diameter of 1,000 X and f/D=2, except 
that P2 has a shorter focal length f/D=l. The scan loss is summarized in Figure 4.19. 
Conclusions are listed below. 

(i) For mechanical scan by tilting reflectors, the best system is PI. The scan loss 
at 0 O =8° (114 beamwidths) is only 0.6 dB (Figs. 4.5, 4.6). The sidelobe 
level for the 0 O =8° position is increased considerably (from -31 dB to -13 
dB). This problem, which may be alleviated by using the cluster 
compensation method [25, 27, 29], should be studied. 

(ii) The folded version of PI is the Cassegrain antenna Cl. In the present study, 
the feed is taken to be a single open-ended circular waveguide with a=3X. As 
a consequence, the amount of folding achieved is small (the length reduction 
is from 2,000 X to 1,650 X). If more folding is desired, a much larger feed 

should be used. 

(iii) To shorten the f/D from 2 to 1, reflector P2 must rely on a cluster feed to 
reduce its scan loss. The excitation of the cluster varies as the beam scans. 
The scan performance of P2 is still not as good as that of PI, indicating that a 
19-element cluster cannot totally compensate for the reduction in f/D. 

(iv) Tilting the subreflector of a Cassegrain antenna can only achieve a small scan 
(about ±15 BW). It can be used in conjunction with the electrically more 
effective but mechanically more costly main reflector tilting to achieve a small 
local scan. 

(v) Among the six antennas, only the offset parabolic reflector P3 scans the beam 
electronically. The price is high since (a) there arc two identical antennas, one 
to scan up and one to scan down, (b) the feed has 65,000 elements, and (c) 



with a 19-element feed cluster, the scan loss at 00=8° is 6.3 dB. Without the 
cluster, the loss is 15.4 dB. This is much worse than the 0.6 dB loss for PI. 



a. Pl : Tilting Parabolic Reflector 



b. P2= Same as PI Except f/D= I and Cluster Feed 



c. P3 : Off-set Parabolic Reflector With Electronic Scan 



|*566X*j 


65 k 

Elements Total 


2a=2.l3X 
Cluster Feed With 
1,7, or 19 Elements 
"On" at Any Time. 



Figure 4.1. Single reflector antenna systems. 





a. 


C I s Tilting Main Reflector 




C3= Tilting Both Reflectors 



Figure 4.2. Cassegrain dual reflector antenna systems. 
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Beam vs. Antenna Efficiency 



Beam efficiency and antenna efficiency for PI symmetrical parabolic reflector 
as a function of feed raduis. The feed is an open-ended waveguide. 


Sym. Parabolic Reflector Untilted 



Theta (deg) 


Figure 4.4. PI symmetric parabolic reflector far-field pattern for 0° scan. 
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Primary Feed Pattern 



. Primary feed pattern for PI symmetrical parabolic reflector. The feed is an 
open-ended circular waveguide. 


Sym. Parabolic Reflector Tilted 4 deg 



Theta (deg) 


Figure 4.6. PI symmetrical parabolic reflector far-field pattern for 8° scan (114 
beam width scan). 
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Sym. Parabolic Reflector Tilted 10 deg 
70 
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Figure 4.7. PI symmetrical parabolic reflector far-field pattern for 20° scan (286 
beam width scan). 
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1 .0 > Magnitude > 0.8 
0.8 > Magnitude > 0.6 
0.6 > Magnitude > 0.4 
0.4 > Magnitude > 0.2 
0.2 > Magnitude > 0.0 


Figure 4.8. Relative excitations for the 1 9-element feed cluster in P2 symmetrical 

parabolic reflector: a. beam scanned 0°, b. beam scanned 8°, and c. beam 
scanned 20°. 
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igurc 4.9. P2 symmetric parabolic reflector far-field pattern for 0° 


f/D=l Sym. Para. Reflector Tilted 4 deg 



scan. 


Figure 4.10. P2 symmetric parabolic reflector far-field pattern for 8° scan. 



f/D=l Sym. Para. Reflector Tilted 10 deg 
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Figure 4.1 1. P2 symmetric parabolic reflector far-field pattern for 20° 




a. 


b. 




1.0 > Magnitude > 0.8 
0.8 > Magnitude > 0.6 
0.6 > Magnitude > 0.4 
0.4 > Magnitude > 0.2 
0.2 > Magnitude > 0.0 


Figure 4. 12. Relative feed excitations for the P3 off-set parabolic reflector, 
a. Beam unscanned. b. Beam scanned 8°. 





Off-set Para. Reflector with 19 El. Feed 



Theta (deg) 


Figure 4. 13. P3 offset parabolic reflector far-field pattern for 0° scan. 


Off-set Para. Reflector with Cluster Feed 



Figure 4.14. Directivity vs. scan for the offset parabolic reflector. Note that as the number 
of cluster elements increases, so does the directivity. 
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Off-set Para. Reflector with 19 El. Feed 



figure 4.15. P3 offset parabolic reflector far-field pattern for 8° scan. 


Cassegrain Reflector Untilted 
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Figure 4. 16. Symmetrical cassegrainian reflector far-field pattern for unscanned beam. 
Cl, C2, and C3 have identical patterns for this case. 
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Casseg. with Main Refl. Tilted 5° 



Cl symmetrical cassegrainian reflector far-ficld pattern with main reflector 
tilted to produce 9.96° scan. 



Theta (deg) 


Figure 4. 18. C2 symmetric cassegrainian reflector far-field pattern with subrcflector tilted 
to produce 1.75° scan. 
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Figure 4. 19. Scan loss of the six reflector antennas shown in Figures 4. 1 and 4.2. 
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5. COMPUTATION AND OPTIMIZATION OF BEAM EFFICIENCY 

When a reflector ant enna is used for a radiometer application, the important figure 
of merit is its beam efficiency, not its directivity. In contrast to directivity, there are 
virtually no published results on beam efficiency for a reflector antenna. Typical beam 
efficiency curves exist only for idealized distribution over a circular aperture [31-33]. 
These are useful for general trend prediction, but do not describe the dependence of the 
exact reflector geometry and its feed. 

In this chapter, several topics concerning computation and design aspects of beam 

efficiency are addressed, namely, 

(i) a discussion of three definitions of beam efficiency, 

(ii) the numerical accuracy problem arising from pattern computation codes based 
on physical optics theory, and a method for overcoming this problem, 

(iii) typical beam efficiency data for parabolic reflectors with a circular horn feed, 

(iv) optimization of beam efficiency by a cluster of feeds. 

5.1. Three Definitions 

Beam efficiency is a measure of how well an antenna manages to transmit its power 
within a prescribed narrow cone. It is defined as the percent of total radiated power 
contained in die main beam 

n = £S2£xlOO% (51) 

Meed 

where Pfeed is the power radiated from the feed of the reflector (not including the feed 
mismatch loss). The unsettled question then is how the main beam is defined. Traditionally 
the main beam is taken as an angular region within a circular cone with a half cone angle 
©beam [31], [18], where 
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8bum = polar angle measured from the main beam (where 0=0°) to the (5.2) 
first null of the pattern. 

The definition in (5.2) is sometimes difficult to apply. For a measured antenna pattern, the 
first null is generally not symmetrically located in all cuts, and may not be clearly 
identifiable. Hence, a more practical definition for the main beam is [18] 

Obeam = 2.5 x 63^3 (5.3a) 

where 

03 <jb = half beam width defined by 3 dB down points. (5.3b) 

The second definition also has a drawback. While optimizing the beam efficiency of a 
reflector antenna, it encourages excessive edge taper so that the 3dB beamwidth is very 
wide. That leads to an antenna design with high beam efficiency but poor directivity and 
poor aperture efficiency. This drawback is due to the fact that 03dB is a moving target, 
which prevents a fair comparison between two patterns with different edge tapers. 

In light of the above discussion, we propose a third definition, namely, the actual 
half beamwidth 03 <jb in (5.3) is replaced by a nominal half beamwidth, which is a 
prescribed number, not a moving target For a given reflector diameter D, the nominal 
beamwidth depends on how much one is willing to sacrifice the aperture efficiency for the 
sake of improving beam efficiency. A good choice that conforms to current practice is 


5 360 

03dB " (DA) 


(5.4) 


This choice was arrived at by choosing an ideal aperture distribution with 20 dB edge taper. 
This distribution (5.4) gives the correct location of the 3 dB down point Therefore, 
definitions 2 and 3 are the same for a reflector antenna having a 20 dB aperture taper. This 
edge taper is chosen because it is a good compromise between low directivity loss (0.9 dB) 
and high beam efficiency (98%, not counting spillover). If a reflector antenna does not 
have a 20 dB edge taper, then definitions 2 and 3 can yield very different results. 



Definition 3, however, is better since it does not lead to the use of patterns with excessively 
wide beamwidths as definition 2 does. 
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52 . The Problem with Direct Computation 

Given the secondary pattern of the reflector, a straightforward method to evaluate 
power in the main beam, denoted by Pbeam* is to integrate the field over the main beam, 
namely, 

p ^=i^JL |E|2da <5 - 5) 

where AQ is the half-cone solid angle. This works if the secondary pattern is known 
accurately. To obtain a 2% accuracy in efficiency, the pattern near the main beam region 
must be accurate within 0.086 dB. Such an accuracy is not attainable by standard pattern 
computation methods, which are invariably based on the physical optics (PO) theory. The 
reason can be traced to the fact that the PO does not conserve power. The finiteness of the 
energy condition is not obeyed by the scattered fields which result from the assumed 
surface fields [34]. Technically, this means that the divergence theorem can not be used in 
deriving integral representations for these fields. If we carried out die integral in (5.5) over 
the entire 4ji space, the power in the pattern generally does not match the power incident on 
the reflector, which is feed radiated power minus the spillover. This is due to truncation of 
the surface currents at the reflector edge. We shall use an example to illustrate our point. 

Consider a center fed parabolic reflector with 
Reflector diameter = 2a = 100 X 
F/D = 2.0 

Feed: TE1 1 mode of a circular waveguide with diameter 1.8 X 
located at the focal point. 

Edge taper = 10 dB 
Spillover loss = 16% (-0.75 dB) 
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Its secondary pattern calculated by a PO code is shown in Fig. 5.1. We integrate this 
pattern to calculate the power within a cone with half cone angle Ocone- The result is shown 
in fig. 5.2. Interpolating Fig. 5.2 conservatively out to 0 = 90° results in a total of less 
than 96% power when spillover loss is added in. This clearly indicates that the secondary 
pattern is not entirely accurate and, therefore, power is not conserved. Unfortunately, there 
does not seem to be a simple fix to make PO-based pattern codes more accurate (adding 
fringe current does not solve the present accuracy problem). Since the error is caused by 
truncation of surface currents at the reflector edge, the error is less for cases in which the 
edge taper is higher. For the same reflector with a 2.4 X radius feed (30 dB edge taper), 

the sum of the radiated power and spillover power is 99%. 

For the pattern in Fig. 5.1, let us compare the three definitions for the main beam: 


0.86° 

using definition in (5.2) 

0.85° 

" (5.3a) 

. 0.90° 

" (5.4) 


which gives beam efficiencies, respectively, of 75.5%, 75.5%, and 75.6%. In general, the 
difference in computed beam efficiency will not vary significantly (i.e., the difference in 
efficiency is only a few percent) for values of O^am calculated using the various 
definitions, since the difference between the two solid cone angles is a thin ring located in 
the neighborhood of a pattern null. When a much larger feed (3.0 X) is used, the respective 
beam efficiencies from the three definitions are 90%, 90%, and 88%. 

5.3. Indirect Computation of Beam Efficiency 

In order to achieve beam efficiency results more accurate than those in Section 5.2, 
an indirect method of computation will be described here. The calculated secondary pattern 
is calibrated by an ideal aperture distribution, whose far-field pattern and power enclosed in 
the far-field pattern can both be obtained via closed-form solutions. If the calculated 
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pattern's power enclosed graph (e.g.. Fig. 5.2) is matched carefully to an ideal pattern, 
then the resulting beam efficiency calculation should be more accurate than for the direct 
computation method. 

The following definition is used for calculating beam efficiency by the indirect 
method. The beam efficiency t] is defined as 

'n := ^lpat^l«ccllfeed (5.6) 

where the various elements are defined as follows: 

ripa,: The calculated secondary pattern is examined for relative accuracy. Causes 
of inaccuracies could include overly sparse FFT grids or other errors caused by numerical 
integration. At this point, inaccuracy due to nonconservation of power by the PO method 
is ignored. It is only really necessary for the pattern to be accurate over the first few 
sidelobes since that is where almost all of the power in the secondary pattern is located. 
The cutoff point for the accurate region is denoted by 0**. Then the pattern efficiency is 
defined as 


P(e bewn ) 

" PO.*) 


(5.7) 


where P(0) is the power inside the cone half-angle 0 (Fig. 5.3). 
ti.cc: The exact definition of ti.cc is 


PC0.cc) 

fi = ' 

Power in the actual pattern 


(5.8) 


Because of the nonpower conservation of the PO-based pattern code, it can be difficult to 
calculate the power in the actual secondary pattern. Here, we propose an approximation. 
An ideal closed-form solution for the far-field pattern is used. The following aperture 
amplitude distribution is assumed: 

0(p).C + (l-O[l-(p/a)2)P 


(5.9) 
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where the parameter C is the aperture edge taper in decibels (ET = 201ogC) and the 
parameter P is used to change the shape of the taper (Fig. 5.4) [4]. A similar distribution is 
used in Equation 2 of [31]. If the feed blockage is negligible, then the far-field integration 
over the aperture produces 

T = 7ta 2 [CA 1 (kasin0) + (1 - QAp +1 (kasin0)] (5.10) 


where 


a„(£) * 2 n r(Pfri)-—— 


(5.11) 


The total amount of power in this far-field pattern can be easily determined. Note that the 
pattern is plotted versus kasin0. For a given antenna, the reflector radius a is fixed and the 
pattern truncates at kasin0 = ka [4]. This can then be converted into the cone half-angle 


theta. Then, the r\tcc defined in (5.8) is approximated by 

Tlacc * T|«cc(ideal) (5.12) 

i"| feed : The feed efficiency, Tjfeed, is a measure of how much of the power radiated 


from the feed hits the reflectors). This term takes into account spillover loss. In some 
previous works it is referred to as feed efficiency factor [31], which is a function of the 
feed pattern and the angle subtended by the reflector dish [35]. For a multiple reflector 
antenna system, would appear once for each reflector. 

As an example of the indirect method described in (5.6), consider the following 
case of a symmetric parabolic reflector with the following parameters: 

Reflector diameter = 2a = 1 ,000 X 

Feed: TE1 1 mode of a circular waveguide with diameter 1 .8 X located at the 


focal point 
Edge taper = 10 dB 
Spillover loss = 18% (-0.85 dB) 
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Its secondary pattern is shown in Figure 5.5. From the calculated pattern, we find ©beam 
0.090° (using (5.4)) and we choose 0** = 4.2°. This leads to a pattern efficiency = 
0.956 (Fig. 5.6). Using the ideal pattern with p = 1.0 and a 10 dB edge taper leads to TW 
= 0.997. The feed efficiency is tlfeed = 0.822. The resulting beam efficiency is t \ = 78%. 
On the other hand, a direct (conventional) calculation using (5.1) leads to tj = 76%. When 
a 2.2 X radius feed is used, the corresponding efficiencies are 88% and 86%, respectively. 
Therefore, as expected, the PO error leads to the efficiency being 1-2% low. A limitation 
of the indirect method is that it can only be used in symmetric situations, since the ideal 
pattern is symmetric. 

5.4. Optimization of Beam Efficiency for Cluster Feeds 
5.4.1. Use of cluster feeds 

In some reflector applications a cluster feed consisting of multiple elements is used 
instead of a single feed dement This cluster feed can be used in two different ways: 

• One-to one Excitation in which only one feed in the array is used to produce 

one beam. 

• Cluster Exdtation in which a group of feeds (normally 7 to 50) is used to 
produce one beam. 

The cluster excitation method can produce a more desirable secondary pattern due to its 
higher number of degrees of freedom. The complexity of the feed is also higher for the 
cluster excitation method. The complexity is outweighed if sufficient improvement in the 
results is obtained. 

By property setting the excitations of the cluster elements, it is possible to optimize 
various performance parameters of the reflector antenna. In the past, this has been done for 
optimizing the directivity for wide-angle scans [17,24,36], compensating reflector surface 
distortion [17], reducing sidelobes [24,29,30], etc. In this section we apply the cluster 
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excitation to optimize the beam efficiency. This optimization is solved by means of a 
matrix formulation similar to that employed in [36]. 

5.4.2. Matrix formulation 

Consider a reflector antenna with N feed elements (Fig. 5.7). When element m is 
excited with a unit amplitude (J m = 1) while all outer elements are idle (J n * 0 for n * m), 
the secondary pattern from the reflector is called the active element secondary pattern of 
element m and is denoted by E m (r). In particular, we are interested in the reference 
polarization of En,(r), namely, 

— jkr 

E ra (r)-R*= E^(e,<|>) , kr->~ (5.13) 

I 

Here, R and C are two unitary vectors describing the reference and cross polarizations. If 

all of the feed elements are excited with weighted amplitudes { Jm}* the composite 

secondary pattern is given by 
N 

QflW-IvtfW (5 - 14) 

m=l 

The radiated power within the main beam of pattern Q(0,<}>) is 

p “- = iJL QQ * d£2 (515> 

m=l n=l 

where Zq = 120rt ohms, dQ = sin0 d0 d<J> is the differential solid angle, and the 
integration is over the main beam. The square matrix B has elements B mn 

where 

B ™=-rfL B'WnsWfda <5.i6) 

•'beam 

The power radiated by the feed cluster may be calculated by integrating over the primary 
partem of the feed cluster or by integrating the secondary pattern over all space. These two 
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integrals are roughly but not exactly equal [35]. For the purpose of this thesis the integral 
over the secondary pattern at infinity, which is more valid, is used. This power is given by 
N N 

P^-XXOA^aVIj (5.17) 

m=l n=l 

where 

E?(8,t>[Ef(e,«0]'d£2 (5.18) 

* *spacc 

The ratio of (5.14) to (5.16) is 


Beameff. = 



(J*) T BJ 

<J*) t aj 


(5.19) 


It is now desired to find the J that optimizes the beam efficiency. 


5.4.3. Optimization of beam efficiency 

Given Hermitian matrices A and B, we wish to maximize 'll in (5.19). It is well- 
known that A can be diagonalized in the fashion 

U 4 AU = D (5.20) 


where 

U = a unitary matrix formed by Uj ,U 2 U N ], 

On.Un = an eigenvalue and its corresponding eigenvector of A, 
D = diagonal matrix with diagonal elements d^ = a„ . 


Substitute 

A = (C) T C, 


i 

= =7 =-i 

where C = D U 


(5.21) 


where use is made of the fact that 
the unknown J, 



Introduce a new unknown K in terms of 
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K = CJ = D 2 U _1 J. 
Then (5.19) can be written as 


( 5 . 22 ) 


K* T FK 

k* t k 


(5.23) 


where 


J. _i 

= = 2 =-l = = 2 

F=D U BU 


(5.24) 


Equation (5.22) is the Raleigh quotient in the matrix theory [37]. The maximum value of ti 


q = Xmax, (5.25a) 

which occurs when 

K ■ V max . _ (5.25b) 

Here, X m« is the maximum eigenvalue of F, and Vmax is the corresponding eigenvector. 


Rom die relation 

F V max = Xmax Vmax 
and (5.23), it can be shown that 

A B W max = XmaxWmax 

where 


(5.26) 


(5.27) 


W m .x=UD 2 V m „ 


(5.28) 


Next the solution of J that maximizes q is desired. From (5.21), (5.24), and (5.27), we 


obtain 


j=uu* t w imx =w imx 


(5.29) 


In conclusion, given a feed element cluster with a known geometry, the maximum 
achievable beam efficiency q is obtained by using the element excitations (within a scaling 
constant) con tain ed in J in (5.28), which depends on the element pattern matrices in 
(5.16) and (5.18). 



5.5. Results and Discussion 

5.5.1. Power as a function of sidelobes 

It is of interest to know what percentage of power is located inside a certain number 
of s idel obes. For example, if 99% of the far-field power is within the first three sidelobes 
for a given far-field pattern, then the integral over the secondary pattern to determine this 
power need only cover the solid angle covering the first three or four sidelobes. For the 
i rifa l pattern described in Section 5.3, the integral is closed form (5.9) and the enclosed 
power can be readily found. The power inside the first, second, and third nulls as a 
function of the aperture edge taper is shown in Figure 5.8. This indicates that for patterns 
with hi gh edge tapers, the total radiated far-field power can be approximated by integrating 
over the cone containing the main beam and first three sidelobes. 

5.5.2. Efficiency of reflector antennas 

The results in Figure 5.8 for the ideal aperture distribution do not show the 
dependence on antenna geometry and feed pattern. In order to study the efficiency of a 
reflector antenna, the actual pattern must be used instead of the idealized pattern in (5.10). 
It is interesting to observe that, for a reflector antenna, high efficiency is achieved by severe 
"underillumination." By underillumination, we mean that both the main beam and part of 
the first sidelobe of the feed pattern arc intercepted by the main reflector. This point will be 
explained further in Section 5.5.3. 

5.5.3. Dependence of beam efficiency on feed size 

Consider a parabolic reflector fed by a single circular waveguide with the following 


parameters: 



Reflector diameter = D = 1,000 X 
f/D ratio: f/D = 1.0 and f/D = 0.4 

Feed: TE1 1 mode of a circular waveguide with diameter 2a at die focal 
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point 

We consider the cases: 

Case 1: Symmetric reflector 

Case 2: Offset reflector (lower edge of reflector is on die parabolic axis) 

The variation of beam efficiency with respect to feed radius is shown in Figure 5.9a for the 
symmetric case and in Figure 5.9b for the offset case. The corresponding results for 
directivity are shown in Figures 5.10a and 5.10b, and for edge taper in Figures 5.11a and 
5.11b. It is interesting to note that 

1. Optimum efficiency and optimum directivity do not coincide. 

2. Near optimum efficiency occurs for edge tapers £ 20 dB but optimum efficiency 
occurs when the reflector is slightly "underilluminated" (Le., the first null 

hits die reflector near its edge). 

Optimum directivity occurs for a 12-15 dB edge taper. 

5.5.4. Scanned beam 

A single symmetric reflector with a single circular waveguide feed is used. The f/D 
ratio is 2. 

Reflector diameter = D = 1 ,000 X 

Feed: TE1 1 mode of a circular waveguide located at the focal point 
Two feed radii are examined. The information for them is given in Table 5. 1 . 


Table 5.1. Information for scanned reflector antennas 


Radius (X) 

Edge taper (dB) 

©beam (2) 

®bcara (4) 

Spillover loss (dB) 

1.8 

10 

0.35° 

0.36° 

-0.85 

2.2 

20 

0.36° 

0.36° 

-0.44 
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Scanning is accomplished by tilting the reflector. When the reflector is tilted by a 0 , the 
beam is scanned by angle 0^ - 2a (Fig. 5.12). This is known as the mirror effect and it 
more than doubles the maximum scanning ability of an antenna when compared to that of 
die usual method of shifting the feed and keeping the reflector fixed. The results are shown 
in Figure 5.13. 

An f/D ratio of 2 is very large but large f/D ratios contribute to much better scanning 
results. As an example of this, a comparison of beam efficiencies is made between the 
2.2 X case above and a third case. This case uses a 1,000 X reflector with f/D = 0.4. A 
feed radius = 0.7 X is used to maximize the on-axis beam efficiency. The nominal half- 
cone angle is 0be«n - 0.36° (5.4) and the first null is at 0be*n = 0.44° (5.2). The spillover 
loss is very low (-0.16 dB) and results in a very high on-axis beam efficiency of 94.7%. 
However, even for small scans, the antenna with the higher f/D ratio has a much higher 
beam efficiency (Fig. 5.14). 

5.5.5. Cluster feed 

A single symmetric reflector is used. The f/D ratio is 0.4. 

Reflector diameter = 2a = 1,000 X 

Feed: Seven-element hexagonal cluster of circular waveguides excited by 
the TE1 1 mode. The feed radius is 0.7 X. Cluster is centered at the 
focal point Element spacing is 1.42 X. 

Spillover loss = -0. 16 dB (on-axis with only die center element lit) 

The excitations of the feed elements are chosen as above to maximize beam efficiency. In 
Figure 5.15, results for scanning by an optimized seven-element hexagonal cluster feed 
(ring of six elements surrounding one element) are compared to those for an antenna with a 
single element of the same type as those used in the cluster feed. The feed element size is 
chosen so that beam efficiency is optimized for a single element in die unscanned case. For 
this reason, at small scan angles the optimal result comes from using very small excitations 



on the outer elements; as a result, the improvement is small. The improvement due to using 
the more complex cluster feed is more evident for wider scans. In this case, the efficiency 
improved by over 10% for some scan angles. 



Directivity (dB) 
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Figure 5.1. Secondary pattern for reflector antenna with f/D =2.0 and a 10 dB edge taper. 

The feed is an open-ended circular waveguide of radius 1.8 X excited by the 
TE1 1 mode. 




Power (% of total power in beam) 



Figure 5.2. Power enclosed in the secondary pattern as a function of cone-angle. The 
discrepancy between the enclosed power and the spillover loss is the PO 
error. 






Field strength (dB) 
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Figure 5.4. Aperture distribution for the ideal feed pattern (5.8) with a 20 dB edge taper 
(C = 0.1) and various values of p. 




Directivity (dB) 



Figure 5.5. Secondary pattern with the region of accuracy, bounded by 0»cc, shown. 


Normalized power enclosed in cone 
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Half-angle theta (deg) 


Figure 5.6. Power enclosed, as a function of the cone angle. Almost all of the power in 
the far-field pattern is in the region of accuracy. 
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(Q 0 » $0) 



Figure 5.7. A reflector antenna with a feed cluster. The secondary pattern due to the 
excitation of element 2 (J2 = 1 and all other Jn = 0) is E2. 


Percentage of power enclosed by nulls 
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Figure 5.8. Power enclosed within the main beam and first two sidelobes for the ideal 

aperture distribution given in (5. 10). The power is shown as a function of the 
edge taper of the feed aperture distribution. Note that die curves level out for 
edge tapers greater than 20 dB. 
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Feed radius (a /X) 


Figure 5.9a. Variation of beam efficiency (using definition in (5.4)) with feed radius for a 
symmetrical antenna with various ftDs. The feed is an open-ended circular 
waveguide. 
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Figure 5.9b. 



Feed radius (a/X) 


Variation of beam efficiency (using definition in (5.4)) with feed radius for 
an offset antenna with various f/Ds. The feed is an open-ended circular 
waveguide. 




Feed radius (a/ X ) 


Figure 5. 10a. Variation of directivity with feed radius for a symmetrical antenna with 
various f/Ds. The feed is an open-ended circular waveguide. 





Directivity (dB) 


92 


Figure 5.10b. 



Feed radius (a IX) 


Variation of directivity with feed radius for an offset antenna with various 
f/Ds. The feed is an open-ended circular waveguide. 




Edge taper (dB) 


o 



0.5 1.0 1.5 2.0 

Feed radius (a/ A. ) 


Figure 5.11a. Variation of edge taper with feed radius for a symmetric antenna with 

various f/Ds. The feed is an open-ended circular waveguide (cutoff occurs 

if the radius is roughly less than 0.4 A). 




Edge taper (dB) 


o 


Figure 5.11b. 



Feed radius (a/ X ) 


Variation of edge taper with feed radius for a symmetric antenna with 
various f/Ds. The feed is an open-ended circular waveguide (cutoff occurs 
if the radius is roughly less than 0.4 X). 
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6 


scan 


* 2a 


Figure 5. 12. The beam is scanned by tilting the reflector and keeping the feed fixed. 



Beam efficiency (%) 
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Scan angle theta (deg) 


Figure 5. 13. Comparison of scanning ability for two antennas with different sized feeds 
(single open-ended waveguide excited by TE1 1 mode). 




Beam efficiency (%) 
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Scan angle theta (deg) 


Figure 5. 14. Comparison of scan-loss results for antennas with different f/D ratios. A 

large f/D ratio allows much better scanning. The feed is a single open-ended 
circular waveguide excited by the TE1 1 mode. 




Beam efficiency (%) 
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Figure 5.15. Comparison showing the improvement obtainable by use of a cluster feed 
with feed excitations optimized for beam efficiency. The antennas have f/D 
= 0.4 and all feed elements are open- ended waveguides with radius = 0.7 A. 
The cluster feed used seven elements in a hexagonal arrangement (six around 
one) with center-to-center spacing = 1.42 A. 




6. USE OF FREQUENCY SELECTIVE SURFACES IN 
REFLECTOR ANTENNA DESIGN 


Frequency selective surfaces (FSS) are very valuable for the design of multiband 
reflector antennas. In a standard antenna design, it is desired to place the phase center of 
the feed at the focal point of the main reflector. However, in multiband applications it is 
often necessary to use more than one feed. In this case the feeds must be kept physically 
separate. This can be achieved by placing one or more of the feeds at an unage point of the 
focal point This is often done in dual reflector configurations, such as cassegraiman or 
gregorian. In a cassegrainian configuration, the subreflector is located between the main 
reflector and its focal point (Fig. 6.1). An additional feed could be placed at the focal point 
of the main reflector, but its energy would be blocked by the subreflector. This problem is 
avoided if the subreflector is transparent to the energy emitted behind it, but reflects the 
energy emitted from the cassegrainian feed. This effect can be achieved with use of an 
FSS. In this chapter, the use of FSS in the design of the ATDRSS triband reflector 
antenna is examined. 

6.1. The ATDRSS Project 

NASA’s Tracking and Data Relay Satellite System (TDRSS) presently provides a 
vital link in space communications. The TDRSS satellites substantially increase earth-to- 
space link availability and provide a near continuous exchange of information. A single 
TDRSS satellite can transmit and receive high-data-rate information to and from low earth 
orbiting spacecraft via two single access (SA) reflector antennas. These steerable SA 
antennas can provide simultaneous S-band and Ku-band communications with one 
spacecraft at a time. Communications to and from orbiting spacecraft can also be 
accomplished via an S-band multiple access phased array antenna, though at much lower 
data rates. A separate space-to-ground link antenna operating at Ku-band provides 



communications between a TDRSS satellite and the TDRSS White Sands Ground Terminal 
in New Mexico. 
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In order to provide additional bandwidth for increased communications demand, the 
advanced TDRSS, or ATDRSS, project has been proposed and is scheduled for launch in 
1997. As conceived, the ATDRSS satellites will incorporate Ka-band capability in the SA 
reflector an tennas, in addition to the S-band and Ku-band services. Therefore, in order to 
pwt these future requirements, the development of a triband reflector antenna for ATDRSS 
is critical. 

6.2. Design 

There are two approaches to designing a triband reflector antenna. Both approaches 
use a single-band feed for S-band and isolate the S-band from the other two bands by 
tramns of a frequency selective surface (FSS). In the first approach a multiband feed 
(usually a corrugated horn) is used. In the second approach two single-band feeds are used 
which are isolated from each other by means of a second FSS. Both approaches have 
pluses and minuses. The multiband feed is more compact and avoids additional FSS losses 
[38-40]. However, design of a multiband feed is much more difficult than that of a single- 
band feed. It is more difficult to optimize parameters such as feed size and taper. 
Performance will further degrade if the phase centers of the two bands do not coincide. In 
contrast, the horn design for the second approach is much easier and feed losses will be 
lower. The design will not be as compact since there is an extra feed and FSS. The 
challenge is in designing a low-loss FSS. The task is made more challenging by the 
relatively small separation (about 2: 1) between the Ka-band and Ku-band. FSSs are today 
typically used to discriminate between bands with a 6:1 ratio (for example between S-band 
and Ku-band). 

In this thesis, we have chosen the second approach, in which three separate feeds 
arc used. It is felt that the improvement in feed performance will not be offset by the higher 
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FSS losses. It is also easier to analyze the feed system in the second approach since the 
design of the multiband feed is to some extent a hardware design problem. Each feed is 
optimized for a single band and the feeds are isolated by means of frequency selective 
surfaces (FSS). The FSS will transmit certain frequencies while reflecting others. Two 
reflector antenna configurations are presented below, an offset single reflector (Figs. 6.2, 
6.3) and a symmetric shaped dual reflector (Figs. 6.4, 6.5). The advantages of each 
design are presented in Table 6.1. 


Table 6.1. Design summary 


Main reflector 

FSS 

Advantages 

parabolic, offset 
Fig. 6.1 

2 planar 

planar FSS only 

smaller diameter main reflector (12.5 ft) 
solid reflector, not mesh 

shaped, symmetrical 
Fig. 6.3 

1 planar and 

l 

1 curved 

smaller volume (shorter focal length) 
similar to existing TDRSS design 


6.2.1. FSS design 

Planar FSSs have an advantage in that they can be theoretically analyzed by 
methods such as Floquet modes. Curved FSSs must be analyzed as being locally flat The 
initial design is then tested to observe the perturbation caused by the curvature. Several 
iterations are then usually necessary to completely compensate for the effects of the 
curvature. 

FSS designs have been developed that will provide the necessary transmission and 
reflection characteristics. The FSSs use ring elements, due to the circular polarization of 
the radiated field. A total of four FSSs were designed. FSS1 and FSS2 are used in the 




offset configuration. FSS3 and FSS4 are used in the symmetric configuration. A 
break do wn of the requirements for each surface is shown in Table 6.2. 


Table 6.2. FSS requirements 


FSS 

S-band 

Ku-band 

Ka-band 

1 

transmit 

— 

reflect 

2 

transmit 

reflect 

transmit 

3 

transmit 

reflect 

reflect 

4 

— 

transmit 

reflect 


FSS3 is the only curved surface, being the subreflector of the dual cassegrainian design. 
Using these designs, far-field patterns have been computed for the above antenna systems, 
including losses due to the FSS effects. 

6.2.2. Offset reflector 

The offset configuration uses a reflector with a diameter of 150", a focal length of 
130", and an offset height of 15". Both FSSs are planar. FSS1 is farther from the 
reflector and is orientated vertically (Fig. 6.3). This surface uses Arlon DICLAD880 that is 
15 mil thick and has a dielectric constant e = 2.17 - j0.0017. The lattice angle is 60° (as it is 
for all four FSSs); the dimensions for the ring elements are shown in Figure 6.4. FSS2 is 
tilted 5° from vertical The substrate is Arlon DICLAD880 but the thickness is 30 miL The 
dimensions of the ring elements are shown in Figure 6.5. All feeds are assumed to be 
corrugated circular horns. The S-band feed has D — 10" (1.8 X at 2.2 GHz), the Ku-band 
has D = 1.6" (2.0 X at 14.9 GHz), and the Ka-band feed has D = 0.94" (2.0 X at 25.25 

GHz). 
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6.2.3. Symmetric reflector 

The symmetric configuration uses a cassegrainian subreflector. The main reflector 
has D = 168" and a hole at the vertex 28" in diameter. The subreflector has D = 28". The 
subreflector is shaped to avoid sending energy into the hole in the main reflector. The main 
reflector is shaped to avoid energy blockage by the subreflector (Fig. 6.6). This shaping is 
subtle and at S-band frequencies the main reflector and subreflector appear to be parabolic 
and hyperbolic in shape respectively. The Ku-band and Ka-band feed structures are 
inside the shadow cast by the subreflector on the main reflector so that feed blockage losses 
will be minimal (Fig. 6.7). The S-band feed is located at the focal point of the main 
reflector. At S-band the main reflector can be considered to be parabolic because the 
deviation due to shaping is only a fraction of a wavelength. This feed is a crossed dipole 
wi th 4.4" diameter subreflector. The Ku-band and Ka-band feeds arc corrugated circular 
horns that are 3.6" and 1.88" in diameter, respectively. These feeds are significantly larger 
than the corresponding feeds for the symmetrical case. The offset design has f/D * 0.87. 
For the dual reflector design, the distance from the feeds to the subreflector is about 1.5 
times the diameter of the subreflector. Therefore, the Ku-band and Ka-band feeds for the 
<hi$ii design need to be more directive to maintain spillover losses comparable to those for 
the offset design. The subreflector is FSS3, transmitting at S-band and reflecting at the 
higher frequency bands. The substrate is Arlon DICLAD810 with dielectric constant £ = 
10.5 - j0.0158 and thickness 14 mil. The dimensions of the periodic element are shown in 
Figure 6.8. FSS4 is tilted 35° from vertical. The subtrate is Arlon DICLAD880 with a 
thickness of 200 mil. There are two layers with the rings facing inward towards each 
other. The separation between the layers is 84 mil and is also filled with DICLAD880. 
The dimensions of the ring elements are shown in Figure 6.9. 

6.2.4. Comparison of reflector size 

Most Phase A designs for ATDRSS have proposed 16* diameter mesh main 
reflectors [41,42], The large diameter compensates for blockage, feed losses, and RMS 
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errors which are higher for a mesh surface than for a solid surface. By combining high 
feed efficiency and a shaped reflector to lower blockage losses, the size of die reflector is 
reduced to 14' for our symmetric design. For the offset case, there is no blockage. An 
assumption of a solid surface reduces RMS surface losses. In this case a diameter of 12.5 
is achievable while meeting link budget requirements. If the solid reflector is hinged to 
allow folding, then it is possible to fit a solid reflector of this size on the launch vehicle 
(Space Shutde or Atlas-Centaur). 

6.3. Results 

A physical optics-based model for an open-ended circular waveguide is used as the 
feed. Physical optics is used to calculate the field incident on the reflector (or main reflector 
in the case of the dual reflector design). An FFT is then used to calculate the far-field 
pattern. 

It is important to integrate the effects of the FSS into die reflector analysis, because 
the transmission and reflection coefficients of the FSS are a function of the incident angle 
of the radiation. The theta and phi components of the incident field interact differently with 
the FSS. Therefore, the incident wave must be broken down into its theta and phi 
components as defined in the FSS coordinate system. The local z-axis is chosen as being 
the normal to the surface at the point of intersection, and directed into the same half-plane 
as die global z-axis (Fig. 6.10). The reflected and transmitted fields are then related to the 
incident field via a matrix formulation. 


H? 


ti «2 

H? 

. Hf . 


. t 3 u. 

.HfJ 

H? ] 


n r 2 

Hf 

. Hf . 


. *3 r 4. 

Hf _ 


and 


( 6 . 1 ) 



In the case of a transparent surface for transmission and a perfect electrical conductor (pec) 
surface for reflection, (6.1) reduces to 
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In the reflector system the energy spreads out as it leaves the feed. Therefore, it is 
incident on the FSS over a wide range of incident angles (Fig. 6.11). For example, in the 
offset configuration the energy from the Ka-band feed is incident on FSS1 at a range of 
incident angles of 6.6° < 0 < 64.8°. The strength of the field also varies as a function of the 
angle from the feed axis. In addition, the losses due to phase shifting by the FSS must also 
be considered. These factors can not be adequately accounted for unless the FSS effects 
are integrated into the reflector system. In general, the FSS losses are less than 1 dB and 
are usually on the order of a few tenths of a dB. The most noticeable effect is reduction of 
the null in the cross-pol at boresight However, in all cases the cross-pol is at least 20 dB 
below the ref-pol. The results for the two configurations are shown in Table 6.3. The 
FSS losses are shown in parentheses. TRW’s estimated link budget requirements [43] are 
also shown for comparison. 


Table 6.3. Directivity results 


Band 

Freq. (GHz) 

Offset dir. (dB) 

Symmetric dir. (dB) 

Required (TRW est.) (dB) 

S 

2.2 

37.2 

(0.1) 

37.7 

(0.1) 

36.0 

Ku 

13.7 

— 


54.6 

(0.6) 

51.0 

Ku 

14.9 

54.0 

(0.1) 

55.4 

(0.5) 

51.0 

Ka 

25.25 

58.7 

(0.1) 

60.2 

(0.4) 

54.0 

Ka 

27.5 

58.8 

(0.6) 

60.3 

(1.0) 

54.0 
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The 2.2 GHz, 14.9 GHz, and 25.25 GHz frequencies are the center frequencies 
(the lower edge of the receive band) for the S, Ku and Ka bands, respectively. Results 
were computed for two additional frequencies where FSS losses were a maximum. For 
both designs, at the upper edge of the Ka band (27.5 GHz), increased FSS losses almost 
completely negated the increase in directivity resulting from using the same-size reflector at 
a higher frequency. For the symmetric design, FSS losses in the Ku band were highest at 
the lower band edge (13.7 GHz). In all other cases, the FSS losses at the band edges were 
lower than or roughly equal to losses at the center frequency. At every frequency, the 
offset design has lower directivity than the symmetric design, despite having lower FSS 
losses. This is due to the fact that the offset design uses a smaller reflector. 

It should be noted that Table 6.3 does not give the complete picture. The computed 
directivities shown include spillover/illumination losses, blockage losses, and FSS losses. 
They do not include reflector surface losses, radome losses, feed losses, and feed line run 
losses. These losses are taken from published Phase A results [41-43]. When these 
additional losses are added, all link budget requirements are still satisfied. The gain margin 
(over TRW’s Phase A report specifications [43]) is shown in Table 6.4. The margin is 


Table 6.4. Antenna gain margin over TRW specifications 


Band 

Freq. (GHz) 

Gain Margin (dB) 
Offset Symmetric 

S 

2.2 

0.7 

0.8 

Ku 

13.7 

2.3 

2.3 

Ka 

27.5 

3.3 

2.9 


computed in each band for the frequency at which the margin is a minimum. This is 2.2 
GHz for S-band and 27.5 GHz for the Ka-band. The actual Ka-band margin is actually 




slightly lower, since the TRW Ka-band specification is for 25.25 GHz, not 27.5 GHz. 
For the Ku-band, the margin is computed at 14.9 GHz for the offset design and at 13.7 
GHz for the symmetric design. 
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6.4. Physical Layout 

The offset fed single reflector antenna configuration has been selected for 
integration into proposed ATDRSS designs. The offset design uses a solid reflector and 
planar FSSs, which can be accurately modeled by existing computer programs using modal 
analysis. The offset fed single access antenna geometry shown in Figure 6.2 has been 
adapted for a conceptual spacecraft design. Some of the assumptions for this design are (a) 
an Atlas Centaur launch vehicle, (b) an Advanced Communications Technology Satellite 
(ACTS) size spacecraft body and ACTS type solar arrays and (c) the S-band amplifier can 
be located in the spacecraft body and the Ku- and Ka-band equipment can be located in the 
antenna arms near the feed horns. 

CAD AM drawings have been completed for the design concept and a 1/1 3th scale 
model has been constructed, employing this offset-fed antenna configuration. Figure 6.12 
is a photograph of the spacecraft model in the launch ready state. Figure 6.13 shows the 
spacecraft model with both offset fed antennas fully deployed and pointed north and south. 
This configuration allows for a full 360° offset antenna scanning capability. Lunar access, 
which has been recently added to the ATDRSS mission, is achievable with this 
configuration. 



Main reflector 



Figure 6. 1 . Cassegrain antenna. The focal point of the parabolic main reflector coincides 
with a focal point of the hyperbolic subreflector. 
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Figure 6.2. Offset design for proposed ATDRSS triband reflector antenna. The single 
reflector dish is parabolic and has a solid surface. The reflector is 150” in 
diameter with an offset height of 15” and a focal length of 130”. 
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Figure 6.6. Dual symmetric design. The reflector is shaped and has a diameter=168”. 
The subreflector has diameter=28". 
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Figure 6.7. Close-up of the symmetric design feed system. The S-band feed is a crossed 
dipole with reflector, and the other feeds are corrugated circular horns. 
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Figure 6.9. Geometry for FSS4. 
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MEW 


may be curved. 



Figure 6. 1 1 . Since the feed docs not emit a plane wave, the energy radiated is incident 
a surface over a broad range of incident angles. 
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NASA 

C-91 


1/13 Scale Model of Advanced Tracking and Data Relay 
Satellite (ATDRS) Concept Stowed for Launch 


2.2m K~band 
downlink antenna 
stowed for launch 


/r 3.8m S, K» and K u -band single access 
I Vantennas folded and stowed for launch 



y— S-band multl- 
/ access phased 
array antenna 



Figure 6.12. Model of proposed ATDRSS satellite showing the launch-ready 
configuration (photo courtesy of NASA). 
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0 - 91 - 0*633 

1/13 Scale Model of Advanced Tracking and Data Relay 
Satellite (ATDRS) Concept Configured for Lunar and 
Out of Ecliptic Plane Missions 



pointed toward lunar 
orbiting satellite 

C041 4*001 


Figure 6.13. Model of proposed ATDRSS satellite showing the triband reflector antennas 
in the fully deployed mode (photo courtesy of NASA). 


7. CONCLUSIONS 


We have expanded existing methods of calculating far-field patterns for reflector 
antennas to include many of the difficulties presently encountered in reflector antenna 
design and analysis. By using methods of analysis that arc generalized and allow for more 
variation, the scope of problem types that may be tackled is broadened. The purpose is to 
develop methods of analysis flexible enough to handle tomorrow’s problems as well as 
today’s. Using these techniques, problems such as spillover loss for reflector antenna 
waveguide feeds, optimization of beam efficiency for reflector antennas, and the analysis of 
reflector antenna systems including frequency selective surfaces have been addressed. 
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Abstract 

A frequency selective surface (FSS) with concentric metal ring elements are analyzed by 
a mode matching technique. The currents on the metal rings are represented by modal functions 
of a coaxial waveguide. Both the circumferential and radial field variations are accounted for. 
As a result, the present analysis applies to thin as well as thick rings. 

Several FSS configurations are designed for the ATDRSS application, which requires the 
separation of three frequency bands: S, Ku, and Ka bands. A typical configuration consists of a 
thin layer of periodic concentric metal rings printed on a 10 to 20 mil substrate. The substrate is 
supported by a three-layer honeycomb structure for mechanical strength. All five layers are used 
in the present analysis. Our calculated results are in good agreement with measurements. 
Conclusions about using FSS for the ATDRSS and similar projects are given below. 

(i) Among various FSS elements (dipole, tripole, cross, rectangular or circular aperture), 
the circular ring element stands out as the best in terms of control of pass-stop bands, and 
stability with respect to incident angle/polarization. Another advantage of rings is that rings of 
different radii can be concentrically arranged for multi-band operations. 

(ii) A theory and a computer code 'cring.f have been developed under the present project. 
The code can be used to analyze, design, and optimize ring FSS. A dozen of good FSS designs 
generated by cring.f are given. Some of them have been verified experimentally. Those designs 
can serve as templates for future applications. 

(iii) FSS designed for Ka-band operation requires very stringent tolerances, namely, the 
geometrical dimension be accurate within 1 mil, and the dielectric constant be precise and 
homogeneous within 5%. If those tolerances are not met, FSS performance deterioted rapidly. 

(iv) Because of the tolerance problem, FSS in the Ka-band can be built with confidence 
only in a planar form, but not in a doubly curved form that is required for a symmetrical 
reflector antenna. As a result, the use of FSS in an ATDRSS-type reflector a is restricted to off- 
set configurations. Of course, the tolerance problem will be eased as mechanical and material 
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controls improve. Future use of FSS in a ATDRSS type system, even for a curved sub-reflector, 
is promising. 

This report consists of three parts. The first two are view graphics used in the project 
review meetings in January and October, 1992. The last one contains the mathematical analysis 
of the FSS with thick ring width. 
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Introduction 

Analysis of frequency selective surface (FSS) is rapidly becoming a mature subject in the 
electromagnetic theory [1-2]. Invariably, FSS is modeled by an infinitely large periodic 
structure. The reflected and transmitted fields in space are represented by discrete Floquet 
modes. The analysis problem then becomes the determination of the Floquet modal coefficients 
by matrix equations. 

There are many different metal elements can be used and have been used to form a periodic 
FSS. Examples are rectangular plate, circular plate, cross, Jerusalem cross, dipole, tri-pole, 
rectangular rings, and circular rings. In the past years, the circular ring element has emerged as 
the best element for applications because of its 

* good pass band and good band characteristics, 

* relatively stable performance with respect to incident angles and polarization, and 

* concentric ring configuration that allows multiple band operation. 

In the last year project, we used the ring elements to design a set of four FSS for the ATDRSS 
application [3]. That design was done with a computer code named 'ring.f . 

In ring.f computer code, there are two assumptions on the structure of the FSS. First, the 
width of the metal rings are assumed thin in terms of electrical wavelength. As a result, the 
current on the ring has no radial component and no radial variation [4]. The second assumption 
is that the metal ring sheet is rested on a single dielectric substrate. 

Both assumptions turn out too restrictive for the ATDRSS application. To order to better 
control FSS performance, it is necessary to use rings with wide width, say 0.4 wavelength width. 
Furthermore, in actual reflector application, FSS is glued to a supporting honeycomb structure. 
It has been found by a study a NASA Lewis that support structure has profound effect on the 
FSS performance, and it must be taken into consideration in the analysis. 

A new computer code named 'cring.f is developed in the current project. The current on 
the metal rings are represented by the coaxial waveguide modes, which have both circumferential 



156 


as well as radial variation and field components. In cring.f, there is no restriction on the number 
of substrate layers. Thus it can be used to model the composite structure of FSS and its 
honeycomb support. A description of the analysis that cring.f is based on is presented in this 
report. New designs for ATDRSS application are given in [5]. 

Reference 

[1] C. C. Chen, "Scattering by a two-dimensional periodic array of conducting plates", IEEE 
Trans. Antennas Propagat., Vol. AP-18, pp. 660-665, 1970. 

[2] S. W. Lee, " Scattering by dielectric loaded screen", IEEE Trans. Antennas Propagat., 
Vol. AP-19, pp.656-665, 1971. 

[3] E. A. Parker and J. C. Vardaxoglou, " Influence of single and multiple-layer dielectric 
substrates on the band spacing available from a concentric ring frequency selective 
surface". International Journal of Electronics, Vol.61, pp.291-297, 1986. 

[4] M.L.Zimmerman and S.W.Lee, "Integrated reflector antenna design and analysis", EM 
Lab. Rept MZ1-91, University of Illinois, Urbana, IL 61801, July 1991. 

[5] S.W.Lee, M. Christensen and Sean Ni, " Multi-ring FSS analysis using coaxial waveguide 
modes", EM Lab. Rept GF92-2, University of Illinois, Urbana, IL 61801, October 5 1992 ( 
this report is included in this volume). 
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Chapter 1. Mathematical Analysis of FSS 

In this chapter, the formulations for calculating the transmission and reflection from an 
infinitely large periodic frequency selective surface (FSS) are derived. The FSS is made of 
capacitive multiple metal rings. The computer code used in the later chapters is based on the 
formulations presented in this chapter. 


1.1. Field Representation 

The geometry of the rings is shown in Figure 1.1. In the mathematical analysis, let us 
assume that the metal ring is not perfectly conducting, but with a finite surface resistivity R c 
(e.g., Rj. = 100 ohm per square). In the current version of the code, Rc is set to zero. As will be 
shown in the analysis, the addition of finite Rc is a relatively minor step. We can add the finite 
Rc capability to “cring” with a small effort. 

The configuration in Figure 1.1 is horizontally stratified. A typical region defined by zj < z 
< Z 2 , the field is derived from two potentials 


TE: E = - V x (z\|/). (1.1.1) 


9w 

1 9 2 \|/ 

H * “ jkoZofi 5x5z 


9\ i/ 

1 9 2 \|/ 

Hy " jkoZoi 1 dyd 2 



e -° 

TM: H=^V x (zip) 

c _ J_ d 2 V |i _J_dv 

x "jkoe^ Hx= Zo 3 y 
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SIDE VIEW 


Figure 1.1. A resistive FSS illuminated by an incident plane wave. 




MODE MATCHING FOR PERFECT CONDUCTING SCREEN 



Figure 1.2. Mode matching equations for a perfectly conducting screen. 



MODE MATCHING FOR RESISTIVE SCREEN 



Figure 1.3. Mode matching equations for a resistive screen. 
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Due to the periodic nature of the problem, the total fields in all regions, except in the resistive 
layer, can be represented by Floquet space harmonics. Explicitly, they are given by 


Transmitted region 

V = II T pqQpq ( x > y) exp(-j7^ z ) (1.1.3a) 

pq 


Dielectric layers 

V = 12 Qpq (X. y) [ A SJ exp(-j-$ z) + B<"> exp(+jT« z)] < 1 . 1 .3b) 


Incident region 

V = II Qpq (x, y) [l pq exp(-jy^pq } z) + R pq exp(+j^ ) z)] (1.1 .3c) 

P9 


Here the summation indices (p, q) take the values 


p, q = 0, ± 1,±2, ... 

The factor Q pq describes the transverse variation of the (p, q)th space harmonics: 

Qpq (x, y) = exp[-j(u po x + x pq y)] (1.1.5) 

where 

u po = (2pJl/a) + k 0 V M-26 £ 26 sin e 0 cos 4>0 ( 1 • 1 - 6 ) 

Vpq = (-2ptt/a tan Q) + (2qtt/b) + lq) V M-26 e 26 sin % cos <j> 0 (1.1 .7) 

The propagation constants in (2.3) are defined by 


1/2 

Tpq = [ko Mn e n - ( u p 2 o + v m)] * n = 1 , 2 , ..., 10 


(1.1.8a) 
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where the square roots are taken such that 



For given incident field (I p q, Ipq), the unknown modal coefficients {Tpq, A^j, Rpq} are to 

be determined. In addition to the TE modes described in (1.1.3), TM modes exist in the 
various other regions. They are derivable from vji in accordance with (1.1.2). We represent \j/ 
in the same manner as in (1.1.3) except that the constants {Tpq, A^j, .... Rpq} are 
replaced by {Tpq, aJJ, ... , Rpq}. 

1.2. Fields at the Interface 

We use the same field representations given in Section 2, 79-Rept. The tangential fields just 
below the R-sheet are (Appendix C) 

E t (z = 0-) = II j W pq Qpq (X, y) { (Otpq [(P^ + P ( p 3 q } ) Ipq + (p^ + p^) Rpq] 

pq 

+ (^W&pq[(Pw , -pS)Tpq + (P < N-P < N , )Rpq] } 2.1) 

H,(z = 0-) = XX j Wpq Qpq (X. y) { (lp! 1 >/k 0 Z 0M6) Ppq [ f Ppq “ P«) [ pq 

+ <p£ - Pffi) Rpq] + H/Zo) OCp, [(pg} + p®) Tpq 

+ (Ppq+Ppq)fipq] ) d-2.2) 

where 

«pq = (X Vpq - y Upo) = -z x Ppq (1.2.3a) 

A -1 .A A x A A 

Ppq = Wpq (X Upo + y Vpq) = Z X 


(1.2.3b) 
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w pq “ po + v pq 

(1.2.3c) 

0pq = ' (Vpq/Up Q ) 

(1.2.3d) 

p = 0, ± 1,±2 ±P 

(1.2.4a) 

q = 0, ± 1 , ± 2 ± Q 

(1.2.4b) 


From (1.2.4) we note that the total number of Floquet space harmonics used in our computation is 
NF = 2 (2P + 1) (2Q + 1) (1.2.4c) 

The tangential fields just above the R-sheet are 

Ef(z = 0+) = XX j Wpq Qpq (x, y) { Otpq (ppq + Ppq) Tpq 
pq s 

+ (-/,i 1 W l )^ N (pS-p < S>Tp< l ) <‘- 2 - 5 > 

H t (z = 0+) = XX j Wpq Qpq (x, y) { (Tpq^oZoHl) Ppq (Ppq — Ppq) ^pq 

+ C-1/Zq) Otpq (pg + pg) Tpq } (1.2.6) 

Next, let us consider the representation of the unknown surface current on the resistive Xi- 
Note that, regardless of the value of R c , the current J s satisfies the same boundary condition that 


Normal component of J s = 0 


(1.2.7) 
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at the rim of the circular plate. In the special case in which R c = 0 (perfectly conducting plate), 
J s is given by the well-known expression 


where 


J s (z = 0) = ^ z x E ind (z = 0) 


( 1 . 2 . 8 ) 


Ejnd (z = 0) = Tangential electric field in the circular aperture of a 
complementary (inducting) screen when R c = 0 


= XXf(C V +D V ) U V + (C^ 1 +D H ) U H 
LL, l v'-mn + u mnJ u mn T v'-mn T u inn^ ^mn 


mn 


r r ^ 

i mn 


V^c J 


(C V - D V ) V V + 

V'-mn u mn/ T mn T 


f T ^ 

1 mn 


M. 


(C H - d h ) v h 1 

V'-mn v mn J 


(1.2.9) 


V'Hl'-C J 


The representation in (1.2.8) is based on the Babinet principle that 


ZoEjnd (z = 0) — » H cap (z - 0+) or H cap (z - 0-) 
and the fact that 

J s (z = 0) = £ x [H cap (0+) - H cap (0-)] 


Because of (1.2.7), the same representation (1.2.8) for the conducting case can still be used for 
the present resistive case. 

The summation in (1.2.9) represents the superposition of fields on four rings. Let us 
concentrate on the field on a typical ring. As explained in Appendix A, the assumption of 
narrow ring width allows us to use a simplified modal field representation, namely. 


[Am “I 1 ' 

2 sin m<(> 

[_ln(d/c)J 

.-cos m<J)_ 


V VH = 0 
v m u 
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In another code “ring2c,” we used the exact coaxial waveguide modes to represent the field on 
the rings. For narrow rings, there is virtually no difference in their solutions. 

Our problem at hand is to determine NF transmission coefficients { Tpq, T pq } and NG 
reflection coefficients { R^, R^ } , for a given set of incidence coefficients { I pq , Ipq } . In the 
present formulation, we determine aperture-field coefficients { C mn + D mn , C mn - D mn } as an 
intermediate step. 

13. Field Matching at z = 0 

There are three boundary conditions to be satisfied. The first one is 

Eton (z = 0-) = E (z = 0+), inZj+12 (1.3.1) 

where 

Xj = resistive circular plate 

X 2 = remaining portion in a unit cell. 


We multiple both sides of (1 .3. 1) by the operator 


a/2 b/2 

f dx r dy exp jCupoX + v pq y) 

-a/2 -6/2 


We obtain 


T = 

A pq 


Ppq 


Ppq 


(6) R pq + 
+ Ppq 


> (,) + d (3) 

J pq Ppq T 

.(5) A „(6) A pq 

>pq 


Ppq 



0 (2) -p (4) 

Ppq Ppq p 

,p(5)_p(6) K M 

LPpq rpq 


. Ppq - Ppq T 

pS-pS M 



which relate the transmission coefficients to the reflection coefficients. 


(1.3.2) 


(1.3.3) 


(1.3.4) 


The second boundary condition is 
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J s (outside) = 


J s (0) , in Xj 

0 , in X 2 


(1.3.5) 


Here J s (0) is the induced current on the resistive circular plate and is given in (1.2.7). The 
current J s (outside) is that calculated from the discontinuity of the tangential H-field outside the 
resistive plate, namely, 

J s (outside) = z x [H t (0+) - H t (0-)], in Xi + X 2 0 -3.6) 

We will now consider the enforcement of (1.3.5) in detail. Substituting (1.2.2) and (1.2.6) into 
(1.3.6) gives 


J s (outside) — XXjWpq Qpq (x, y) { otpq Fpq [(Ppq + Ppq) ^pq (Ppq + Ppq) ^pq] 


pgi v , ♦ L *8 KpS - pg) ♦ <Pg - pg) fiJ + ^ ) 


where 


p(3) _ 

r pq 


-iy^)pg-pg,ifTg>g-pg 

2 JfkoH.jpg) + - 2 + pW 


4 [*£-*£] 
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,(2) rt (4) 


pq 


Ppq-I 


2 w 


pq 


X 


(i) 

pq 





p£ + p (4) 


pq 


5 (2) -d (4) 

Ppq Ppq J 


1_J_ 
2 W 


pq 



(1.3.9) 


(1.3.10) 


(1.3.11) 


In deriving (1.3.7), we have made use of (1.3.3) and (1.3.4). We multiply both sides of (1.3.5) by 


We obtain 


(1.3.12) 


a/2 b/2 

" A “ 

a pq 

I dx I dy [expj(UpoX + v pq y)] [Wpq] 

A 

a/2 — tv2 

-Ppq - 


ab Ww { [(P^ + p$> !pq + (Pm + Ppq) R pq] ^M + Ww ) 


- yy T x ( c ^ + d v ^x (3)H + d** ) 

— L~ ^pqmn v^mn ^ ^pqmn VMnn T ^mn/ 


mn 


_v ( r _ n V fC H - )1 

^pqmn V^mn ^mn/ ^pqmn v^mn u mn/ J 


(1.3.13) 
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j ab { [(p® - p®) Ipq + (p® - p ( *>) Rp,] F® + I„F® } 


-tt rx (I)v ir v +d v ux (1|B rc H + d h i 1 

— ***** L yv pqmn V'-'mn T u mn/ T -' v pqmn V'-'mn T ^mn/ J 


where 


f d)V 

v pqmn 

jf(l)H 

>-' v pqmn-i 


Y ( 3 > V ‘ 
yv pqmn 

Y (3)H 

‘- yv pqmn- 1 


-iw CT (,) 

J w pq u pqmn 


-iw a (3) 

J w pq u pqmn 


pqmn 

cos m (Jfpq 
Lsinm^pq. 

sin m 4>pq 
L-cos m <J> pq . 


V (4)V 
^pqmn 

yP 

i-^pqmn-J 


= -J w pq 


/fe 'N 
mn 


V^c J 


t (4) 

■'pqmn 


sin m <t>pq 
|_cos m ((Jpq. 


J (1) 
J pqmn 

J (3) 

■'pqmn 


given in Appendix B 


j (4) =0 

J pqmn u 


In deriving (1.3.13) and (1.3.14), we have made use of the fact that 
z x A *«p q = -A • p pq 
z x A ‘Ppq = A • apq 

— V V 

From (1.3.13) and (1.3.14) we can solve for {Rpq, Rpq) in terms of fCi^p, 


(1.3.14) 


(1.3.15a) 

(1.3.15b) 

,...}. The results are 
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Rpq- 


Ppq + Ppq T F pq T 
'pg + pffi ^ 


, ( ^pq) YY Jv(3)V /pV |-vV \.y(4)V /pV _ qV \ 

+ ? / (?] mv “ 1 ^pqmn v^mn ^ u mn/ ^ ^pqmn v'-mn mn/ 

jabw^(p® + pW) mn 


X (3)H (f .H H \ . x (4)H /pH _pH , 1 
^ ^pqmn v^mn T u mn/ x ^pqmn V^mn ^mn/ J 


ftOJ-nC 3 ) p(4) 

- _ Ppq Kpq - r Pq T 

M "pg-5S "'sgflsg-PSS) " 


+ : — = XX { X^ v „ (Cl + d!) + x$£ (Cl + Dl) } 
jab W p q (p pq P pq ) mn 


So far we have expressed the transmission and reflection coefficients in terms of unknown 


current coefficients 


r pV pjV pH pjH pV j-jV ^-,H pH | 

l ^mn» u mn’ '-mn’ ^mn' ^mn> u mn> '“mn’ u mn J 


(1.3.18) 


Those expressions are independent of the resistivity R c . 

The third and the last boundary condition to be enforced is 

E t (z = 0-) = R c J s (0), inlj (1.3.19) 

where E t is given in (1.2.1), and J s in (1.2.8). Multiply both sides of (1.3.19) by 

A A » tV 

d 2ji • z X U mn 

J pdp f d<J) (1.3.20) 

o o .g xV v - H 
L ^ A ▼ mn -I 
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we obtain 


XX i 

pq 


-x 


(3)V 

pq/w 


, -V(3) H , 

I- A pq/w-J 


Rpq + 


v k 0 


rd)V 

^pq/w 


rim 


L-^pq/w-l 


R pq 


= (-j)(-D M r)Z 0 


C /w + D /w 


r H 4 . n H 

■Ww + u /w- 


(1.3.21) 


where 


pq 


r 

r y (4)v i 


< 

-^pq/w 

R ' ^ 
R pq 


Y (4)H 


• A -pq/w-l 

> 


r-Hx- d m tiz 0 + 


1 /w 
V Mc y 


*~/w ^/w 


r H _n H 

.Ww ^/w. 


R W = (ppq + pg) Ipq + (Ppq + pg) R P q 
R w = (pg - pg) Tpq + G>g - Ppq) Rpq 


(1.3.22) 


Substituting (1.3.16) and (1.3.17) into (1.3.21) and (1.3.22), we obtain the desired matrix 
equation. 

%’C' = G' (1.3.23) 

Here 


X' = 


l 

<— ► 

t (l)H 

Y<i)V 

•jKDH " 

t (3)V 

t (3)H 

j<3)V 

t<3)H 

t (5)V 

t (5)H 

j(5)V 

Y<5)H 

_ t (7)V 

t (7)H 

f 7)V 

t<7)H _ 
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^mn + ^mn 


"£ r 

c£„ + d“„ 


& 


G' = 


^mn “ ^mn 



-cl - 51. 


-gffi- 


The explicit expressions of X' and G' are given below. 


1 /wmn 


,(3)V ,H 
- /wmn - 





_L Y ( 3 > V ’ H rr( 3 > 

„(3) Pqmn u pqmn 

F pq 


sin m 

cos m <l>p q _ 


+ 



\ y(l)V,H 

( 3 ) A pqmn 

F pq 


o (1) 

cos m <|)pq 

> 

u pqmn 

_sin m <t> pq _ 

j 


+ 67 6^ (-i y-i T|Z 0 


1 /wmn 

= Same as 

- t (ir 

except < 

t( 3)V,H 
_ 1 /wmn J 


_t (3) _ 

< 


t (5)V ,H~ 
1 /wmn 

= Same as 

-fsr 

except < 

t (7)V ,H 
- l / wmn J 


_t <7)_ 




V 


x (3) -+ x (4) 

> 

x (1) ->o 

no diagonal term- 
X (4) -+ X (3) 

> 

no diagonal term-* 


(1.3.24a) 


(1.3.24b) 


(1.3.24c) 
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t{5)V,H‘ 
1 /wmn 


t(7)V,H 
. 1 /wmn . 






(4)V,H (4) 
pqmn °pq/w 


sin m (^pq 
_cos m <t)pq_ 


> 


+ S? Sw (-ly- 1 TiZo(r /w /koe c ) - (1.3.24d) 


■«sr 

= 11 4 

r 

sin m (j> pq 


pq 

< 

cos m <t> pq _ 


w o^, (F^ 4 ^ / F^) I 
w pq W pq/w V 1 pq ' 1 p>q/ A pq 


cos m 4»pq 


Wpq <7, 


pq/w (-F < w'F < m)('>S /k o e 6> I w 


Lsin m <))p q 


g* 

1 ^ 

= XI 4 

r 

sin m <J> pq 

> 

si 

00 

1 

PQ 

s. 

_cos m (Jjpq. 

J 


In summary, the final matrix equation is in (1.3.23) for the unknown coefficients {C'J. Once 
{C'} are solved, the reflection coefficients are calculated from (1.3.16) and (1.3.17), and the 
transmission coefficients from (1.3.3) and (1.3.4). 
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Appendix A: Modes on a Narrow Ring 


To determine the modal field on a metal ring, we may consider its complementary problem, 
the modal field on a coaxial waveguide, because these two fields are related by the Babinet 
principle. Consider a coaxial waveguide shown in Figure A.l. Its normal modes can be found in 
a standard textbook (e.g., N. Marcuvitz, Waveguide Handbook , New York: McGraw-Hill, 1949, 
pp. 72-80). They are given in terms of Bessel functions and are rather complicated. Under the 
condition of narrow gap that 


d - c 
c 


« 1 


(A.l) 


we can use a set of simple and approximated modes to be deduced below. 

Consider the field in the neighborhood of Q in the coaxial waveguide. Under the condition 
(A.l), it can be well approximated by that in a slot waveguide on the right-hand side of Figure 
A.l. The fields in the latter are expressible in terms of two vector potentials, namely. 


TE modes: 


f 2 m 7t ,1 


y = A 


mn 


COS 


2nd 


. 2 m 7t 
L sm 2nd x J 



(A. 2a) 


TM modes: 


V = A mn 


cos 


2 m 7t 
2 7t d 



sin 


2 m 7t 
2 7t d x - 



(A. 2b) 



where 
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The boundary condition for the slot waveguide is that 


(Field at AB) = (Field at A'B') 


(A.3) 


Figure A. 1 . The field in a coaxial waveguide with a narrow gap can be approximated 


by that in a slot waveguide. 
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which leads to the admission of both cosine and sine functions in (A.2). Making use of the 
replacement 

x' — > d <(> (A.4a) 

y' — » p - 0.5 (d + c) (A.4b) 


we obtain the approximated modes in a narrow coaxial waveguide, namely, 


\\f~ A 


mn 


cos m <)> 
[sin m <() 


cos ^ (p - c)| e ± -’ I ' mnZ 


(A. 5a) 


cos m (J) 


• .1 ±jr mn z 


■ sin ~(p - c)) e' 
sinm<t>Jl o 


V“A mn | 

Because of the fact that d -» 0, we make a further approximation by using only 


(A. 5b) 


n = 0 


(A.6) 


modes in (A.5), we obtain 


TE modes: 


V- 


yi 

Au m 

m=0 


cos m <b 
sin m (J) 


[C m e- jr ™ z + D m e +jr ™ z ] 


(A.7a) 


TM modes: V|/ = 0 


(A.7b) 


where 


(A. 7c) 


The non-zero field components can be calculated with the formulas 
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E 


P 


l_cty 

P 


(A.8a) 


H 1 1 d 2 V 

* j ko Zq p d<t> az 


(A.8b) 


Returning to Figure A. 1 let us write down a field representation at z = 0 for the coaxial 
waveguide: 


E,(0) =1 { (Cl + D V J + (c£ + DS) u£} 


(A.9a) 


where 



(A.9b) 


U 


V,H 


m 


A m Zo 

1/2 

sin m 6 

_2rt In (c)_ 


-cos m 4>_ 


,m = 1,2, 3, ... 
,m = 0, 1,2, ... 


(A. 10) 


where Z 0 = 120 It and A m = 1 if m = 0 and A m = 2 if m * 0. The m = 0 mode is the well-known 
TEM mode in a coaxial waveguide. Note that U^’ H has been normalized such that 



U^ H pdpd<t> = 


Zq, if m = n 
0, if m * n 


(A. 11) 


is satisfied. 



177 


Appendix B: Coupling Coefficients 

To change (1.3.5) into a matrix equation, we apply the operation (1.3.12) to both sides of the 
equation. The resultant equation is given in (1.3.13) and (1.3.14). In this appendix, we give the 
details of deriving the coefficients (Xpq mn ). 

Applying (1.3.12) to the right-hand side of (1.3.5), the integral is reduced to 4 integrals over 
the 4 metal rings, namely, 


f -If 

Jover cell axb A * J c 


4 rings 


over a ring 


In the following, let us concentrate on a typical ring with outer and inner radii (d,c): 


r Y (i)v,Hi 

A pqm 

= ff pdpd(| l c i( V l + W>w po 

&pq 

V (3)V,H 

J —slot 

R 

L-^pqm J 


LPpqJ 


■U, 


V,H 


m 


(B.l) 


where U^’ H are given (A. 10), and 

&pq = (x v N - y u po ) = -z x (3 pq 
Ppq = Wpi (x Upo + y Vpq) = Z X tt pq 

The final results after evaluating (B.l) are 


cS 

X 

1 

tl) 

cos m <})pq 

y (1) h 

L' v pqm- 

— J w pq ^pqm 

sin m <t> pq _ 


(B.2a) 


(B.2b) 


(B.3) 



178 


where 


X 


(3)V 
pqm 

(3)H 

L-'^pqm-i 


-iw o (3) 

J w pq w pqm 


sin m <t> pq 
[-cos m <t>pq. 


4> pq = tan- 1 


V U P°> 


(B.4) 


y (1) = 


A m Zq 


2tc In 


-|l/2 


w pq d 


c ) J 


(-2m 7t) j m J 


w pq c 


(B.5) 


a (3) = 

u pqm 


[27t In (c) J 


-|l/2 


(2k) j m Wpq [J m (Wpqd) J m (WpqC)] 


(B.6) 


The derivation of (B.5) is given below. From (B.l) we have 


f (i)Vl 

w pqm 

AW 

L^pqm J 


= cjdpjd<})e j 


w pqP cos(4> - <]> Dq ) 




sin m <{> 
[-cos m <J)J 


•(-1) Wpq sin(<|> - <t>pq) 


pq' 


(B.7) 


where 


C = 


A m Zp 
2 k • (In d - In c) 



Using the identity in (B.7) 


sin m <|) = cos m <]>pq • sin m(<(> - <|)pq) + sin m 4»pq • cos m(<|> - <J) pq ) 
cos m (j) = cos m <t> pq • cos m(4> - <J> pq ) - sin m 4^ • sin m(4> - 4> pq ) 

we obtain 


rv (1)V i 

^pqm 

y (1) h 

L^pqm-I 



cos m <])pq 
sin m <t> pq 


[cos(m-l) <(»' - cos(m+l) <j>'] e 1 w pq p cos ^ 


(B.8) 
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where 

<t>' = <t> - <t>pq 


The integral with respect to <J)' can be written in terms of Bessel functions, namely, 



which is identical to (B.3) and (B.5). 
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Appendix C: Field Expression in (1.2.1) 


The field just below the R-sheet is the field at z = 0- in region 6 (Figure 1.1). The modal 
coefficient in region 6 is { A^j, } . By matching the fields at the dielectric interfaces between 

layers 6 to 26, we obtain 


pq 


A£ = Ppq + P'w R] 

r ( 6) _ n (3) T n (4) R 

D pq - Ppq L pq + Ppq ^pq 

Apq = Ppq Ipq + Ppq F-pq 

5(6) _ pjO) 7 , r(4) R 

D pq “ Ppq A pq + Ppq R pq 


(C.l) 


To define the coefficients { p^, p^, . . . } , let us introduce the 2 x 2 matrices, for n = 1 , 2, . . . , 26, 


F$ (z) = 


• (n) 
e - J Vt z 


• (n) 

e + JW 


L y^e jY w z -y^ e +jY w z 


(C.2a) 


Fpq (Z) = 


■ (n) 
e “ JY pq z 


. (n) 

e +J V) z 




. .00 


L e n 


(C.2b) 


Then the coefficients in (C.2) are given by 

= fl F^ +1) (/ n ) (C.3) 


The same equation holds for { Ppq. Ppq’ — 1 exce P l that F pq is re P ,aced b y F pJ . 


o (1) o (2) 

Ppq Ppq 


»(3) 


P (4) 
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The field in upper region 1 can be expressed in terms of that in upper region 5, namely, 


A pq = Pw T pq “Pl+jVi' Pel (C4a > 

Bm = P'w T pq PJ (C-4b) 

Aplq = Tpq expl+j-/^ 1 T c ] (C.4c) 

Bpq = p ( $ Tpq CXp[-j7pq X c ] (C.4d) 

where 

x c = metal thickness = zero in the present capacitive FSS 

Ppq = Ppq ex P[”j7pq x c\ ( C 5a ) 

pg = P ( p 9 q } cxpl+hg x c ] (C.5b) 

r p 2 pa i 

= n[F^(/ n )]- 1 F^>(/ n ) (C.50 

.(9) .(10) n 


The same equation (C.5) holds for { pp q \ pg, . . . } except that is replaced by F^. Note 
that { pg, p^q\ p^j, p^ } do not enter into the field calculation. The reason that the 
exponential factors are introduced into (C.5), is to facilitate the field matching at z = T c later. 
The matrices in (C.3) and (C.5c) can be written in a more explicit form, namely, 


r , , . (n) 

(l+ag)e JC P9 

[Fg(/n)]" 1 Fw 1) ( / n) = 5 

L (l-ag)e _jb Pq 

/ x *u( n ) 

(1 _ a pq) e +jb P9 
(l+ag)e +jc m _ 

(C.6a) 

where 



a (n) _ Jn+D-fn) 
d pq ~ Tpq ' tpq 


(C.6b) 



(C.6c) 
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(C.6d) 


n = 1, 2, 3, 4, and 6, 7, 8, 9. 


(C.6e) 


To calculate [F ( p n q } (/ n )] 1 F^ n q +1) (/„), only the factor a^ in (C.6a) is replaced by 



apq (e n /e n+i) 


(C.7) 





